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1. INTRODUCTION 




The purpose of the NASA Aircraft Noise Prediction Program (ANOPP) is 
to predict noise from aircraft, accounting for the effects of the aircraft 
characteristics, its engines, its operations, and the atmosphere. The 
approach to this problem has been on a fundamental basis, as depicted in 
figure 1. The aircraft follows an arbitrary flight path in the presence 
of an observer on the ground. During this operation, noise sources on the 
aircraft emit radiation with defined power,, directional, and spectral dis- 
tribution characteristics, all of which may depend on time. This source 
noise propagates through the atmosphere, being attenuated, to the vicinity 
of the observer. The observer receives the noise signal from the direct 
ray plus a signal from a ray reflected by the local ground surface. 

A number of approaches are available for this general prediction 
problem. These approaches are divided in ANOPP into four categories, 
called functional levels, which are depicted by the schematic in figure 2. 
The functional levels are defined by the amount of data which is processed 
and by the degree of approximation in the prediction methods. Level I 
predicts an effective measure of noise which depends on the observer loca- 
tion. Level II predicts a noise level which depends on the observer and 
time. In Level III, frequency effects are predicted in addition to the 
effects of observer and time. In Levels II and III, the noise measures 
may be subdivided as to the noise source which generates them. Level IV 
predicts the same information as Level III but with more detail in the 
spectral data. 

The present documents deal with Level III. Part 1 describes the 
modules which prepare the data which are required for source noise pre- 
diction and the modules which are used to process the predicted noise. 

Part 2 contains the source noise prediction modules. The principal 
modules used in Level III are depicted in the functional diagram of fig- 
ure 3. Three stages of computation are shown in this diagram. 

The input preparation stage has modules for preparing tables of data 
which are later interpolated by the prediction and output processing 
modules. The Atmospheric Module establishes the properties of the 
atmosphere: pressure, density, humidity, and related variables as a 

function of altitude- The Flight Dynamics Module represents a two-degree- 
flight dynamic analysis wherein the power setting variable is 
used as a control function with a definite time dependence. A second con- 
trol function, the angle of attack, is also used in this module. These 
control functions may be taken from standard procedure tables or may be 
the result of an optimization process where many flights are analyzed to 
find the best operational procedure for a given aircraft performance cri- 
terion. The Source Parameters Modules prepare source parameters for the 
noise prediction modules. These modules take the Mach number, power 
setting, and angle-of-attack functions of time from the Flight Dynamics 
Module and interpolate the engine state table to give engine state vari- 
ables as a function of time. These data are then passed to the second 
stage of computation. Data for atmospheric absorption effects are pre- 
pared by the Atmospheric Absorption Module. This table is used in the 
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Propagation Module in the second stage of computation. The Geometry Mod- 
ule takes data on the aircraft position and orientation from the Flight 
Dynamics Module and evaluates the vectors from the source to each ground 
observer as functions of time. Each noise source will be given in a speci- 
fied axis system such as engine 1 axes or aircraft wind axes so that each 
observer vector will be expressed in several source coordinate systems at 
the same time. 

Noise predictions are made in several coordinate systems which are 
attached to the aircraft. The airframe noise predictions are made in the 
wind axis system by the Aircraft Noise Module. Wing noise has a dipole 
directivity in this system with the Z-axis being the dipole axis as shown 
in figure 4. Predictions are made for virtual observers at a fixed dis- 
tance r and prescribed angle 6 in this sytem at a sequence of times t 
which are widely spaced. Using the wind axis system, which has the axis 
of symmetry, allows for minimum use of the prediction modules since the 
virtual observers are located by only the single coordinate 0. After 
making the predictions for wing noise for a few virtual observers and 
times in the wind axis system, the Propagation Module interpolates these 
data, using the rapidly varying true observer vectors given by the 
Geometry Module, to give free-field unattenuated predictions at actual 
observer positions. The Propagation Module then adds air absorption and 
ground effects by interpolating a previously prepared table from the 
Atmospheric Absorption Module to give actual noise at the observers due 
to the wing noise component from the Aircraft Noise Module. Other noise 
source components in the Airframe Noise Module may be added to the wing 
noise at the source if they are conveniently expressed in the wind axis 
system. Otherwise, the looping process depicted in figure 3 is used to 
propagate the noise to the observer before summation takes place. 

The engine axis systems depicted in figure 3 contain noise source 
components, such as jet and fan noise, with a common axis of symmetry, 
the engine axis. Noise sources are conveniently summed at virtual 
observer locations in these systems before they are processed through the 
Propagation Module and taken to actual observers. Figure 5 illustrates 
this process. Within the engine axis system, modules for fan, combustor, 
turbine, and jet noise are called to make predictions at virtual observer 
positions. Each component noise prediction is expressed at the same set 
of virtual observer positions so that direct summation of the noise may 
be made before propagation to the observers. This summation saves exces- 
sive use of the computational procedures in the Propagation Module. 

Output processing, the third stage of prediction process, involves 
operations on predicted noise at the observer. These predictions are 
functions of frequency, time, source coordinate system, and observer. 

The output processing modules reduce these data by taking summations and 
integrals over these four variables. The Noise Levels Module removes the 
frequency dependence by converting spectra into weighted noise levels 
such as A-level and D- level to further reduce the amount of data. Per- 
ceived Noise Level, a nonlinear measure of noise, is also computed by 
this module. Time integrals are made in the Effective Noise Module to 
evaluate variables such as Effective Perceived Noise. 

An ANOPP Level III noise prediction module is characterized by the 
prediction of 1/3-octave-band noise. The band centers are based on 
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observer frequencies and are independent of time. All other inputs to the 
prediction modules are time dependent. The vectors from the source to the 
observer are naturally dependent on time and observer so that the output 
from a source module is a function of frequency , time^ and observer, in 
that order. 

The prediction of 1/3-octave-band noise is a serious limitation which 
should not be passed over quickly. Some of the most important noise 
sources are actually tones, for example, from the fan rotor of a bypass- 
type engine. In the prediction modules, these tones are assigned to a 
1/3-octave band and subsequently treated as broadband noise. This will 
cause errors in the prediction of atmospheric attenuation, ground effects, 
and even noise levels. Nevertheless, the added complexity of carrying a 
separate procedure for tones suggests that this is not an appropriate task 
for ANOPP Level III and this type of analysis has been relegated to the 
Level IV manual. 

Input and output of the functional modules may be a combination of 
dimensional and dimensionless variables. Two systems of dimensions are 
used, the preferred SI system and alternate U.S. Customary system. One 
principal distinction must be made about the use of the U.S. Customary 
system. The unit of pound-mass is not permitted and must be replaced by 
the less popular unit of slugs. The reason for this is that F = ma only 
if the mass is expressed in slugs. Variables involving heat energy are 
used within the program. These variables, such as Cp, the specific heat 
at constant pressure, are always expressed in mechanical units rather than 
calories or BTU. This facilitates the dimensional balancing of equations 
and the formation of dimensionless groups . 

Within each module the essential computations are made by using 
dimensionless groups. This approach eliminates extraneous variables, 
accounts for scaling, permits the use of an empirical data base, and 
facilitates the use of different systems of units. When the dependent 
data have been generated in dimensionless form, they may be converted to 
dimensional variables for output if this is desired. In many cases, it 
is not desirable to convert to dimensional units since this would increase 
the amount of computation and data. For example, dimensional engine data 
are often tabulated as a function of altitude, velocity, and engine speed 
(or some power setting variable) , whereas it is only necessary to tabulate 
dimensionless data as a function of Mach number and power setting. The 
use of dimensionless groups thus reduces the three-dimensional arrays to 
two-dimensional arrays, a great savings in computation and storage of 
data. 


The sections of this manual are organized along the lines of the 
diagram in figure 3. In each section, a module is described by briefly 
stating its purpose and function, its input and output data, and the 
method or computation process. References are given to the appropriate 
technical literature from which the methods are developed. The module 
description, like the module, is written as an independent section so that 
some duplication naturally results. On the other hand, since modules must 
interface, it is necessary to refer from one module to another in the out- 
put sections. All variables which are not naturally dimensionless will 
have dimensions shown with SI units first followed by U.S. Customary units 
in parentheses. Variables which have been converted to dimensionless 
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groups through division by a reference variable are given with their con- 
ventional dimensional description followed by re and the group for the 
reference variable. As an example, the dimensionless group (T/T-^) 
would be described as 

T* absolute temperature, re 

and, in the symbols list, would be found with a description 

a 

absolute ambient temperature, K (^R) . 

Symbol lists are provided with each section to aid in clarifying all 
nomenclature . 
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Figure 1.- Overall prediction methodology. 
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Figure 3.- Continued. 
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Figure 3 Concluded. 
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Figure 5. Engine component noise prediction for a given engine axis. 








2. AIRCRAFT FLIGHT DYNAMICS 




2.1 ATMOSPHERIC MODULE 


INTRODUCTION 

Selection of an atmospheric model is the first step in predicting 
aircraft noise. Atmospheric properties affect the performance of the air 
craft, the noise generated by the aircraft and its engines, and the propa 
gation of this noise through the atmosphere. Since noise is most signifi 
cant during terminal operations (take-off and landing) , this model is 
only representative of the atmosphere below 10 km. 

Four approximations are used in this model. The first approximation 
is that acceleration due to gravity is constant and is assigned the sea 
level value. This assumption is usually made in the analysis of aircraft 
performance. The second approximation is to neglect the effect of water 
vapor on the pressure, density, and viscosity of air. Consequently, the 
molecular weight of air is taken as the standard sea level value (scale: 
q 12 = 12.0000). The third approximation is that geopotential altitude is 
equal to geometric altitude. For each of these assumptions the maximum 
relative error is of the order 10“^ at 10 km. The fourth approximation 
is that atmospheric properties are functions of altitude only. All equa- 
tions and constants used in this model are taken from references 1 and 2. 

The primary purpose of the Atmospheric Module is to generate dimen- 
sionless arrays of ambient atmospheric parameters at equal increments of 
altitude. Dimensionless parameters are used to provide simplified forms 
for the equations and to facilitate working in different systems of units 


SYMBOLS 

c speed of sound, m/s (ft/s) 

g acceleration due to gravity, m/s^ (ft/s^) 

H altitude (geometric assumed equal to geopotential) , m (ft) 

Ah output altitude increment, m (ft) 

h absolute humidity, percent mole fraction (ratio of number of 

H 2 O molecules to total number of molecules in mixture, 
1.178 percent at 70 percent relative humidity and standard 
pressure and temperature) 

k coefficient of thermal conductivity, W/m-K (BTU/f t-s-°R) 

M molecular weight of air 

n number of output values 

p atmospheric pressure. Pa (Ib/ft ) 
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R 


universal gas constant, kg-m^/K-s^ (slug-f t^/^R^s^) 


rh relative humidity, percent 

S Sutherland's constant, 110.4 K (198.72°R) 

T temperature, K (^R) 

y = Mgj.(H - Hi)/rTj. 

Y ratio of specific heats 

y coefficient of viscosity, kg/m-s (Ib/ft-s) 

p density, kg/m^ (slugs/ft^) 

Subscripts : 

i input array values 

j output array values 

r standard sea level value 

1 ground level value 

Superscripts : 

* dimensionless value 

average 


INPUT 

Input for constructing the atmospheric model consists of Ah, the 
desired output altitude increment; and P]_, the ground level altitude 

and pressure; and a table of temperatures and humidities as functions of 
altitude. Input of one atmospheric set at ground level consisting of 
altitude, temperature, and relative humidity results in a constant atmo- 
spheric model. Input of more than one set results in a hydrostatic atmo- 
spheric model generated according to the preset altitude increment AH. 

The input altitudes do not have to be supplied at equal increments; how- 
ever, the ground level altitude must be supplied. Output altitudes are 
generated at equal increments. Table I gives the recommended ranges and 
the default values for the input parameters. 


Input Constants 
ground level altitude, m (ft) 

Ah output altitude increment, m (ft) 

2 

p pressure at ground level. Pa (Ib/ft ) 
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Atmosphere Input Table 


H altitude, m (ft) 

T(H) temperature, K (^R) 

rh(H) relative humidity, percent 


OUTPUT 

The output is a table of dimensionless pressures, densities, tempera- 
tures, humidities, sound speeds, average sound speeds, coefficients of 
viscosity and thermal conductivity, and characteristic acoustic imped- 
ances as a function of altitude. 


y 

P*(y) 
P*(y) 
T*(y) 
h(y) 
c* (y) 
c* (y) 

y*(y) 

k*(y) 

P*c*(y) 


Atmospheric Properties Table 
altitude, Mg^(H - 
pressure, re p^ 
density, re 
temperature, re T^ 
humidity, percent mole fraction 
sound speed, re c^^ 
average sound speed, re c^ 
coefficient of viscosity, re 
coefficient of thermal conductivity, re 
characteristic acoustic impedance, re 


METHOD 

The generation of an atmospheric model requires that numerical values 
be established for certain primary constants relevant to the Earth's 
atmosphere. These primary constants are given in table II. All computa- 
tions within the atmospheric model module are based on dimensionless 
equations formed with these primary constants. The basic equations used 
in the computational sequence are: 

Perfect gas law: 


p = pRT/M 


(la) 



which in dimensionless form becomes 


p* = p*T* (lb) 

Equation for speed of sound: 

(p- - YRT/M (2a) 

which in dimensionless form becomes 

(c*)2 = T* (2b) 

Hydrostatic equation: 

dp/p = -Mg^ dH/RT (3a) 

which in dimensionless form becomes 

dp*/p* = -dy/T* (3b) 

where 

y = Mgj.(H - Hj^)/rTj. (4) 

Figure 1 is a graphic representation of the atmospheric coordinates. 

Using the atmospheric variables supplied as input, along with the 
basic equations, the atmospheric model is computed as follows. An alti- 
tude vector Yj is defined for incremental altitude changes Ah by 

Ay = Mgj. Ah/rTj, (5) 


with 


Yj = (j - 1) Ay (j - 1/ 2, n) (6) 

where n is the number of output altitudes. 

The dimensionless temperature is computed by 

T* = T(H^)/Tj. (7) 

Output values for temperature Tj are then found by linear interpolation 
with respect to y. 

Dimensionless pressures are computed by integrating equation (3b) . 

If temperature is assumed to vary linearly between y^ and ^be 

integration may be carried out exactly to yield the recurrence formula: 
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(8a) 


or 



C * * 

Tj ^ Tj.i. 

j ~ 2 ^ 3 ^ 


. . . , n) 






where 


* 



Pl/Pr 


(8b) 


(8c) 


Once the dimensionless temperature and pressure vectors are prepared, 
the other required dimensionless atmospheric vectors are computed as 
follows : 

Density : 


★ 



* 




Sound speed: 


-k 

c . 


] 




(9) 


(10) 


Coefficient of viscosity: 


^ 1, 38313 

^ T* + 0.38313 
3 


Coefficient of thermal conductivity: 


1. 77385 (Xj 

^ -0.0416 /t* 

T + 0.8516 X 10 '3 

3 


Characteristic impedance: 


* * 


P 



( 11 ) 


( 12 ) 


(13) 


The constant 0.38313 is the ratio of S/T^, where S is Sutherland's 
constant . 
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Average sound speed is defined by 


c(y) = ^ (14) 

r y 

I rc(y)1 ^ dy 

•'0 ^ 

where the denominator of equation (14) is the time for vertical transmis- 
sion- If temperature is again assumed to vary linearly between yj and 

yj+]_f the denominator of equation (14) can be integrated exactly to yield 

the following recurrence formula: 


where 


I . 
3 



+ 


2 (Ay) 



(15a) 



with the following condition: 


Yl = 0 


(15b) 


(15c) 


Equation (14) is then computed in dimensionless form for altitude yj by 


c 


(y .) = y ./I . 
^3 ^3 3 


(16) 


The input humidity is expressed as relative humidity in percent. For 
computational purposes, it is more convenient to express humidity in abso- 
lute terms as the mole ratio (in percent) of H 2 O molecules relative to the 
total number of molecules in a mixture. The absolute humidity is defined 
in terms of temperature, pressure, and relative humidity by 

, , 8.4256-(i0.1995/t*)- 4.922 log T* 

hj = (rhj/pj)l0 ' ' 3-' : (17) 


where rhj is computed by linearly interpolating with respect to y. 

Once the atmospheric values are computed in dimensionless form, 
dimensional values for printed output are computed using the following 
equations : 


Hj = (RTj„/Mgj,)yj + ( 18 ) 

Pj = PrP5 <19) 
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( 20 ) 


T . = 


Pj = PrPj 


Cj = C^Cj 


- 

c . = c c • 
D ^3 


Pj = yj.pj 


kj = k^kj 


PCj = Pj-C^PCj 


( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 


REFERENCES 

1. U.S. Standard Atmosphere, 1976. NOAA, NASA, and U.S. Air Force, 

Oct. 1976. 

2. Sutherland, Louis C.: Review of Experimental Data in Support of a 

Proposed New Method for Computing Atmospheric Absorption Losses. 
DOT-TST-75-87, U.S. Dep. Transp., May 1975. 
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TABLE I.- RANGE AND DEFAULT VALUES OF INPUT PARAMETERS 


Input 

parameter 

Minimum 

Default 

Maximum 

Ah, m 

1 

100 

100 

, N/m^ . . . 

90 000 

101 325 

110 000 

H / , m • • • • 

-300 

0 

10 000 

T, K 

200 

288.15 

300 

rh, % 

0 

70 

100 


TABLE II.- STORED PRIMARY CONSTANTS 


Constant 

SI Units 

U.S. Customary Units 



9.806 65 m/s^ 

32.1741 ft/s^ 

M 

28.9644 

28.9644 

R 

8314.32 m^/K-s^ 

49 718.96 ftV°R-s^ 

Y 

1.40 

1.40 

Pr 

1.013 25 X 10^ Pa 

2116.22 Ib/ft^ 

Pr 

1.225 kg/m^ 

0.002 377 slug/ft^ 

Tr 

288.15 K 

518.67 °R 



340.294 m/s 

1116.45 ft/s 

Pr 

1.7894 X 10-5 kg/m-s 

3.737 X 10~^ slug/ft-s 

kr 

6.0530 X io~6 w/m-K 

4.0674 X io“® BTU/ft-s-°R 

RT^/Mg^ .... 

8.434 515 6 X iq 5 m 

2.767 210 65 x lo'* ft 
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Altitude (H^) 


Altitude 


Ground Level 


Standard 
Sea Level 
(H-0) 



1.- Dimensional and dimensionless altitudes 
for the atmospheric model. 
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2.2 GEOMETRY MODULE 


INTRODUCTION 

The calculation of aircraft noise requires a complete spatial and 
temporal definition of the noise propagation paths. The geometric rela- 
tionship between the noise source components and the observers must be 
defined for the duration of the aircraft flight. The aircraft position is 
given, either from input or from the Flight Dynamics Module (FLI) , in an 
Earth-fixed coordinate system as a function of time. The aircraft orien- 
tation is also given with respect to the Earth-fixed coordinate system in 
terms of Euler angles. These coordinate system conventions are described 
in detail in reference 1 and in other standard texts on aerodynamics of 
the airplane. Although the Flight Dynamics Module uses two-degree-of- 
freedom equations, the Geometry Module will be valid for six degrees of 
freedom so that ANOPP may be extended to the more general case of three- 
dimensional flight paths by input of a three-dimensional path or by addi- 
tion of a more general flight module. 

Observer positions are given as input to the program in an Earth- 
fixed coordinate system. The task of the Geometry Module is to find the 
vectors from the noise source on the aircraft to the observers and to 
express these vectors in terms of a coordinate system at the aircraft. 

These vectors are functions of time because the noise source moves with 
the aircraft. 

It is not desirable to make noise predictions for each observer at 
all times during the flight. An observer close to the point of brake 
release will have large values of noise at different times than the 
observer several miles down the flight track. Noise at an observer which 
is more than 20 dB below the maximum value is of little interest in most 
cases. Consequently, a significant savings in computation results when 
the times for which the noise is certain to be negligible are eliminated 
in advance. 

The noise emitted from an aircraft takes a finite time to reach the 
observer. Since noise predictions are desired at evenly spaced values of 
the reception time, the reception times are treated as an independent 
variable and a set of flight emission times are determined. The aircraft 
position and orientation are then found at these emission times for each 
observer. Intermediate reception time points are selected at sufficiently 
small increments to ensure that the changes in noise characteristics 
between successive time points are not large. 

This module produces a table of the emission time, distance, polar 
directivity angle, azimuthal directivity angle, and elevation angle to 
each observer as a function of reception time and source coordinate system. 
This table is used by the Propagation Module (PRO) for the computation of 
propagation effects. 
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SYMBOLS 




c 

AdB 


h 

-V ^ 
if jfk 

k 

M 




(1^1 , ^ 3 ^ 


s 

T 

T 

-^r 

t 

At 


(x,y,z) 


y 

Y 

9 

A0 

P 




reference area of aircraft, (ft^) 
speed of sound, m/s (ft/s) 

limiting noise level, down from peak level 

2 

standard acceleration due to gravity, m/s 
observer height, m (ft) 
base vectors 
constant 

molecular weight of air 
initial mass of aircraft 
direction cosines 
observer index 

universal gas constant, m^K-s^ (ft^/^R-s^) 
distance, m (ft) 
source coordinate system name 
transformation matrix 

standard sea level temperature, K (^R) 
time , s 

reception time increment, s 
position coordinates, m (ft) 
altitude, re RTj./Mg^ 
elevation angle, deg 
polar directivity angle, deg 
polar directivity angle limit, deg 
density, kg/m^ (slugs/ft^) 
azimuthal directivity angle, deg 
Euler angles, deg 


(ft/s^) 
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Subscripts: 


a aircraft or Earth-fixed flight axes 

i time index 

m local minimum 

min global minimum 

o observer or Earth-fixed observer axes 

r standard sea level 

s source coordinate system 

w wind axis 

Superscripts : 

* dimensionless value 

average 


INPUT 

Input to the Geometry Module consists of the aircraft body axis posi- 
tion and orientation vector. It is understood that these data are given 
with respect to the Earth-fixed flight axis system of figure 1, The 
orientation vector is taken to be the set of three Euler angles described 
in reference 1. The average speed of sound as a function of altitude is 
obtained from the atmospheric properties table. There may be several 
noise source axis systems located within the aircraft body system. The 
orientation of these systems are given as input with respect to the body 
axis system. All source coordinate systems are assumed to have the same 
origin. Observer positions are given in the Earth-fixed observer system 
shown in figure 1. Table I gives the recommended ranges and the default 
values for the input parameters. 


Input Constants 

A^ reference area of aircraft, m^ (ft^) 

AdB limiting noise level, down from peak level 

k constant 

m^ reference mass of aircraft, kg (slugs) 

t^ initial time, s 

final time, s 
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At reception time increment, s 

A9 maximum polar directivity angle limit, deg 


Flight Dynamics Table 


t 


Xa(tr 

ya(t) 

Za(t) 


flight time, s 


aircraft body coordinates, m (ft) 


^{tr 

Oa(t) 


aircraft body Euler angles, deg 


^ (t) 
w 

0w(t) 


aircraft wind axis Euler angles, deg 


Observer Table 


o 


observer index 


Xo(o) 

Yo(o) 

Zo(o)_ 


observer coordinates, m (ft) 


Source Coordinate System Table 


H'3(s) 

0s(s) 

$s(s) 


source system ncune 


source Euler angles, relative to body axis, deg 
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Atmospheric Properties Table 


y altitude, re RT^^/Mg^ 

c (y) average speed of sound, re 


OUTPUT 

Output from the Geometry Module is the vectors from the source to 
the observer. These are expressed as spherical components in the source 
coordinate system. In addition, the elevation angle from the observer to 
the source is provided for use in the Propagation Module. 


Geometry Table 

t^ reception time, s 


o 

observer index 


s 

source coordinate system 

name 

tg(tQ,0,S) 

emission time, s 


r (tQ,o,s) 

distance, m (ft) 


6 (tQ,o,s) 

polar directivity angle, deg 

(J) (tj^,o,s) 

azimuthal directivity 

angle , 

Y (1^,0, s) 

elevation angle, deg 



METHOD 

Each observer has an associated time segment during which the emitted 
noise dominates over other times during the flight. This time segment 
(or segments) is to be estimated based on relative distance to the 
observer . This estimation can be done by scanning the vector from the 
source for each observer as a function of time, recording its minimum 
magnitude, and then eliminating those times for which the magnitude is 
greater than a fixed constant, say 10, times the minimum value. This 
operation will define a limited time slice for each observer and further 
computation of noise for that observer may be made within that time slice 
only. 


Elimination of Excess Data 

First, the flight times from the flight dynamics table are examined 
and all times less than the initial time t^ and greater than the final 
time t^ are discarded. In addition, if the points are closer together 
than the characteristic flight time increment At^ given as 
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( 1 ) 


At 


f 




the extra points are discarded. The constant k in equation (1) is a 
user parcimeter. 


Supplementing Sparse Data 

Input to the Geometry Module may consist of data at widely spaced 
time points. If the time points are spaced at time intervals greater than 
the input interval limit At, the aircraft body coordinates, body Euler 
angles, and wind axis Euler angles are interpolated with a cubic spline to 
produce intermediate values spaced at intervals less than At. 


Location of Minimum Distance 

The computations now proceed separately for each observer, and the 
process continues until all observers have been completed. The observer 
coordinates are transformed into the Earth~fixed flight coordinate system 
by the axis reflection transformation: 


r— — 1 


“l 0 0 



o 




o 

yo 

= 

0-10 


Vo 

_Zo_ 


0 0 -1 


h 


The position of the observer relative to the aircraft is then 


X (t) 


X 



oa 


o 


a 

yoa(t) 

= 


- 

Ya(t) 



_^o_ 




Let the flight times be designatedjDy an index i which ranges from 1 
to I and let the observer vector r(tj^) = r^^ be 






(ti) 


+ jyoa<ti) 




(i = 1, 2, 


I) (4) 
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In addition define a discrete path vector by 


= 


- ^i+1 


(i = 1, 2. I-l) (5) 


These vectors, shown in figure 2 , are used to find a minimum distance to 
the observer and the time at which the minimum occurs . 

To find the minimum distance, the first segment is checked to see if 
the projection of r^^ on Sj^, that is, point P in figure 2(b), lies 
within the length of Sj^. This is expressed by 


0 




( 6 ) 


where n^ and n^ are unit vectors in the direction of r and s. If 
this condition is not satisfied, then there is no minimum distance within 
the segment. The equality is included in order that the case of a mini- 
mum at the end point will be acceptable. For each segment which passes 
the test of equation (6) , a local minimiim and time is computed by 


V,i = ^^s,i ^ (^i ^ (7) 

^,i = ti + - ^i) (8) 


For each local minimum, a time point ^ is added to the aircraft posi- 
tion, aircraft orientation, and observer vector data by linear interpola- 
tion on time . 

Finally, the set of values r^ ^ is examined to select the global 
minimum. Then the global minimum distance time 

are recorded for the observer. In the event that there is no local 
minimum, then there can be no global minimum so that the program issues a 
warning message. 


Determination of Time Segments 
The maximxmi distance of interest is given by 


^max 


AdB/20^ 


min 


( 9 ) 


on the basis that spherical spreading alone will give a relative noise 
reduction of AdB in comparison to the noise at the closest point of 
approach. It is desired to find those points r^ which satisfy the 
condition 
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^min 


< 


r • — r 
1 max 


( 10 ) 


and their associated times. This condition is depicted in figure 3. The 
process is one of examining the magnitudes rj^ to find the starting time 

^i “ ” ^max ^ ^i+1 (10a) 

and the ending time 

r. < r(tf) = (10b) 


as depicted in figure 4. Note that the set of values rj^ must include 
the minimum distance vector and at least two other vectors r as depicted in 
figure 5, in order for both equation (10a) and (10b) to have possible 
solutions . 


In order to find the starting and finishing times, each segment is 
examined in turn to find if it contains a limiting point il. As shown in 
figure 6, the equation for the end point can be written as 


a . s . - r . = r 

111 1 1 max 


( 11 ) 


where 


a . 
1 



(i = 1, 2, I-l) (12) 


and tg may be either a start or an end time. Equation (11) gives 


a . = 
1 


n.r.) ± Jft.r.f - s^fr^ - | 

\ii/ ^\ii/ i\i max/ 


(13) 


For the first segment, the condition is added that a^ < 1. This means, 
as depicted in figure 6(a), that the circle with radius l^^tersects 

the vector S]^ or its extension behind time If a < 0, then the 

starting time is before t^^ as shown by equation (12) . For all inter- 
mediate segments, it is required that 

0 ^ ^ 1 (i = 2, 3, 1-2) (14) 
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as shown in figure 6(b). It is possible that both roots will occur and 
define a complete time segment as shown in figure 6(a); however # it will 
generally turn out that only one root of equation (13) satisfies the 
limits in equation (14) . For the final time, the condition 0 <_^a allows 
the final time to lie on the extension of the last path vector 


Limit on Directivity Variation 

The^directivity angle between the path vector s and the observer 
vector r must be limited to a change A0 within each time step. 

The directivity angle is 


0‘ = Arcsin n . ^ n 
1 r,i s,i 


(15) 


and 


A0. 

X 




0 . 

1 


If A0^ > A0, then intermediate time points must be added. 


(16) 


Output Computation 

The direction cosines from the source to the observer in the Earth- 
fixed flight coordinate system are now computed. The position of the 
observer r was computed by equation (3) . Then the direction cosines 
are given by 


•^r,l 

^r,2 


r/r 


L^r,3j 


(17) 


The average speed of sound is known as a function of altitude y 
from the Atmospheric Properties Table. This variable is converted to be 
a function of the flight time by setting 


y(t) 




(18) 


and 


c (t) 


= Cj,c(y) 


* 


(19) 
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The 


where ^ is the direction cosine computed by equation (17) . 
value of reception time t^ that corresponds to each value of t^ can 
now be computed by 

t^ = t^ + r/c(t) (20) 

and a table of t^ as a function of t^ is created. The reception times 
are naturally grouped into segments by the time segments previously 
defined. The reception times found by this process will naturally occur 
at uneven intervals. 

The elevation angle y from the observer to the source can now be 
computed as 


Y = Arcsin n 

' r,3 


(21) 


where n - is the direction cosine. 

r , j 

For computation of the directivity angles ^ the direction cosines of 
equation (17) must be rotated to the source coordinate system. This 
transformation is accomplished with the use of Euler angles for each 
coordinate system as defined in figure 7. 


"r.l 

''r,2 

= [t(<I>s)][t(03)] [t('1'3)] 

"r,l 

^r,2 


s 

-r,3_ 


where 


( 22 ) 


T(vF) = 


cos ^ 
-sin 4^ 
0 


sin ¥ 0 

cos ¥ 0 

0 1 


(23) 


T(0) = 


cos 0 
0 

sin 0 


0 

1 

0 


-sin 0 
0 

cos 0 


(24) 
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1 


0 


0 


T($) 


0 


0 


cos 0 
-sin $ 


sin $ 
cos $ 


Then, the polar directivity angle 6 is given by 
9(t^/0,s) = Arccos (j^r,l)s 
and the azimuthal directivity angle (j) is given by 


(25) 


(26) 


4>(to,o,s) = Arctan (nr,2/"r,3)s 


(27) 


The preceding procedure is repeated until the outputs for all 
reception times, time segments, observers, and sources have been computed 
and the geometry table is complete. The aircraft body axes and wind axes 
are always used for source coordinate systems. Additional source coordi- 
nate systems may be provided by the user . 


REFERENCE 

1. Etkin, Bernard: Dynamics of Atmospheric Flight. John Wiley & Sons, 

Inc . , c . 1972 . 
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TABLE I.- RANGE AND DEFAULT VALUES OF INPUT PARAMETERS 


Input 

parameter 

Minimum 

— 

Default 

Maximum 



1 

1 

1 X lo'* 

AdB 

10 

20 

30 

k 

0.1 

1 

10 

n>o, kg 

1 

416.8 

1 X iq6 

s 

0 

0 


^ 


1 X 10^ 

1 10^ 

At, S 

0.1 

0.5 

2.0 

A0 , deg 

5 

1 

10 

20 
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Figure 1.- Coordinate system definitions for description of 
source to observer geometry. 


OBSERVER 


(a) Discrete representation of flight path. 



(b) Conditions for minimum distance in path segment 
Figure 2.- Location of minimum distance. 
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Figure 5.- Minimum set of three vectors for determination 
of single time segment. 
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(c) Final segment. 

Figure 6.- Solution for starting and finishing 
times for time segments. 


2.2-18 








2.3 FLIGHT DYNAMICS MODULE 


INTRODUCTION 

Accurate prediction of the noise produced by an aircraft in flight 
requires a detailed knowledge of the aircraft position as a function of 
time. The dynamic behavior of any aircraft is a complex problem. A com- 
plete description of the problem requires the solution of six ordinary 
differential equations for the three force components and the three moment 
components in three-dimensional space. The six unknowns for this problem 
are the three position coordinates and the three attitude coordinates 
(Euler angles) of the aircraft. This theory can be found in typical 
textbooks such as references 1 and 2. 

Several approximations can be made to simplify this problem; however, 
each simplification reduces the generality of the description of the 
flight. The first simplification is to restrict the flight trajectory to 
two dimensions. This reduces the problem to three ordinary differential 
equations, two force equations and one moment equation, eliminating any 
consideration of yaw and roll effects. A second simplification is made by 
assuming that all forces act through the center of mass of the aircraft so 
that all moments about the lateral axis of the aircraft are zero. This 
assumption is valid if the aircraft is constantly maintained in a trimmed 
condition. The problem is now reduced to simultaneously solving two ordi- 
nary differential equations for the unknown longitudinal and vertical 
position of the aircraft as a function of time. These approximations are 
used within ANOPP for computation of the aircraft trajectory. If a dif- 
ferent flight trajectory is required, the data can be inserted directly 
as a function of time. 

The Flight Dynamics Module computes the two-dimensional trajectory of 
an aircraft as a function of time. From the module inputs of control set- 
tings, aerodynamic coefficients, engine performance, atmospheric proper- 
ties, and initial/final conditions, the two force balance equations are 
solved for the aircraft position, attitude, and Mach number as a function 
of time. 


SYMBOLS 



reference area of aircraft-engine 
m2 (ft2) 

inlet-face cross section. 


aircraft wing reference area, m2 

(ft2) 

b 

wing span, m (ft) 



aerodynamic drag coefficient, — 

D 

- 


1 2 
2 
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C 


L 


c 

D 


H 

h 




M 


m 






m. 


T 

T 

r 

t 

V 

W 


(x,y,z) 

a 


landing-^gear drag coefficient 

L 

aerodynamic lift coefficient, 

speed of sound, m/s (ft/s) 
aerodynamic drag, N (lb) 
ground force, N (lb) 

acceleration of gravity, m/s^ (ft/s^) 

Heaviside step function 
local altitude, m (ft) 

absolute humidity, percent mole fraction 

landing-gear position 

aerodynaimic lift, N (lb) 

aircraft Mach number, V/c 

aircraft mass, kg (slugs) 

engine air flow rate, kg/s (slugs/s) 

engine mass flow rate, kg/s (slugs/s) 

engine fuel flow rate, kg/s (slugs/s) 

number of engines 

thrust, N (lb) 

standard sea level temperature, K (^R) 
time, s 

engine specific thrust, m/s (ft/s) 
aircraft velocity, m/s (ft/s) 
aircraft weight, N (lb) 

aircraft position in Earth-fixed coordinates, m (ft) 
angle of attack, deg 
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5f flap control variable, deg 

e engine inclination angle, deg 

0 aircraft Euler angle, deg 

0 inclination of flight path with respect to horizontal, deg 

y dynamic viscosity, kg/m-s (slugs/ft-s) 

7T power setting 

p air density, kg/m^ (slugs/m^) 

T coefficient of rolling friction 

Subscripts : 
a ambient 

b body axis 

e engine 

n final limit 

0 break release 

r standard sea level 

1 initial value 
Superscript: 

* dimensionless value 


INPUT 

The amount of information required to define a flight trajectory is 
extensive because of the complex nature of the problem. Description of 
the aircraft geometry must be provided. Initial and final conditions are 
necessary to solve the differential equations and to terminate the solu- 
tion in a finite time. A control varicdsle table defines the aircraft 
attitude and engine setting. An aerodynamic coefficient table and an 
engine performance table describe the characteristics of the airframe and 
engine. Finally, the atmosphere is described by an atmospheric properties 
table. Default values for the input parameters are given in table I. 


Input Constant 

AV maximum velocity increment for engine variable table, m/s (ft/s) 
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Aircraft Configuration Table 


w 


N 


T 


mi 


yi 

hi 

01 




t 

a(t) 

6f (t) 


2 2 

aircraft wing reference area, m (ft ) 

2 2 

engine inlet reference area, m (ft ) 
wing span, m (ft) 

fully loaded reference mass, kg (slugs) 

number of engines 

altitude at break release, m (ft) 

coefficient of rolling friction 

engine inclination angle for each engine, deg 

Initial/Final Conditions Table 
initial mass, kg (slugs) 
initial time, s 
initial velocity, m/s (ft/s) 
initial distance from origin, m (ft) 
lateral position, m (ft) 
initial altitude, m (ft) 
initial flight-path angle, deg 
final time, s 
final velocity, m/s (ft/s) 
final distance from origin, m (ft) 
final altitude, m (ft) 

Control Variable Table 

time, s 

angle of attack, deg 
flap setting, deg 
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ft) 

ft 

power setting for each engine 


^£g 

landing-gear retraction time, s 



Aerodynamic Coefficient 

Table 

a 

angle of attack, deg 


«f 

flap angle, deg 


hA> 

wing height to span ratio 


c^(a,6f 

,h/b) drag coefficient 


(ex f 6 ^ 

,h/b) lift coefficient 



l) landing-gear drag coefficient 



Engine Performance Table 

71 

power setting 


M 

aircraft Mach niomber 


m* (M,7T) 

air flow rate, re D c A 
•^a a e 


m* (M,7T) 

fuel flow rate, re p c A 

a a e 


t*(M,TT) 

specific thrust, re c 

a 



Atmospheric Properties 

Table 

* 

h 

altitude, ^a^r^/^"^r 


c* (h*) 

speed of sound, re c^ 


P* (h*) 

density, re 


y*(h*) 

dyncunic viscosity, re 


h^(h*) 

absolute humidity, percent mole fraction 


OUTPUT 

The Flight Dynamics Module produces two tables that are used within 
ANOPP. The first is the flight trajectory as a function of time. The 
output times for this table are the ones used in integrating the differ- 
ential equations. This table is expressed as a full six-degree-of- freedom 
array, even though the two-degree-of-f reedom assumption is made, to be 
compatible with the Geometry Module. The second is the engine variable 
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table as a function of source time for use in the source parsuneters 
modules. The output times for this table are the control variable input 
times augmented, if necessary, to adequately define the data. 


Flight Trajectory Table 


t^ flight time, s 


(x (t) ,y|^, z (t) ) aircraft coordinates, m (ft) 
(O,0b(t),O) aircraft body axis Euler angles, deg 


(O,0(t),O) 


aircraft wind axis Euler angles relative to body 
axes, deg 


Engine Variable Table 


M(ts> 

^a ^^s ^ 
^a ^^s^ 


source time , s 

Mach number 

engine power settings 

flap setting, deg 

landing-gear position 

cimbient density, kg/m^ (slugs/ft^) 

ambient speed of sound, m/s (ft/s) 

ambient dynamic viscosity, kg/m-s (slugs/ft-s) 

absolute humidity, percent mole fraction 


METHOD 


The two-degree-of- freedom flight trajectory used by this module is 
governed by two ordinary differential equations. The sum of the forces 
along the flight path must be zero and the sum of the forces normal to 
the flight path must be zero. These equations are expressed as 


and 


N 

Fg + L - mg^ cos ® ^ ^ ^e^ ” ^ 


e=l 


N 


-TFg ^ cos (a + Eg) - mg^ sin 9 - D = m 


dt 


e=l 


( 1 ) 


( 2 ) 
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The definition of the force terms in equations (1) and (2) are shown in 
figure 1 for the aircraft during the ground roll and in figure 2 for the 
aircraft in flight. The term V d9/dt in equation (1) is the centrifugal 
acceleration normal to the flight path and the term dV/dt is the 
acceleration along the flight path. 

The ground force term is nonzero only during the ground roll. 

In addition, the centrifugal acceleration and the flight-path angle 0 
are zero during the ground roll. Therefore, applying equation (1) during 
the ground roll yields 


r 


N 


mg 




r ■ L - ^ Tg sin (a + e^) 


(h = ho) 


e=l 


(3) 


V 


(h > hg) 


where h is the local altitude. The coefficient of friction T is a 
function of the landing-gear and surface characteristics. The surface is 
assumed to be uniform throughout the ground roll. It is further assumed 
that the friction force results solely from the main landing gear so that 
it remains constant during aircraft rotation. Since most of the aircraft 
weight is supported by the main gear, the error introduced is small. 

The total thrust T^ for each engine is related to the specific 
thrust and the mass flow rates by 

• it 

T = m c t (M,TT) ( 4 ) 

e e a e 


where 


m = 
e 


m 


+ m. 


(5) 


The aircraft lift L and drag D is computed from the lift and drag 
coefficients as 


L = i P3v\[cL(a,6f,h/b)] 

and 


(6) 


° = J PaV2A^[Cp(a,6f,hA>) + Cjj^j^g(CL)] 


(7) 


In general, the lift coefficient Cj^ and the drag coefficient Cp are a 
function of angle of attack a, flap setting 6f, and altitude h/b. 
Figures 3 and 4 are examples of the effect of flap setting on Cl and Cp 
and figures 5 and 6 are examples of the ground effect on Cl and C^. 
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There is an additional source of drag which is present when the 

landing gear is extended. Figure 7 demonstrates the relationship between 
the landing-gear drag coefficient and the lift coefficient. 

The mass of the aircraft changes as a function of time due to the 
consumption of fuel by the engines. The rate of change of the aircraft 
mass as a function of time is given by 

- I! “fW-’Te) IS 

e=l 

Finally, the position of the aircraft as a function of time in Cartesian 
coordinates is determined from 


= V cos 9 
dt 


(9) 


and 


dz 

dt 


-V cos 0 


( 10 ) 


The numerical solution of differential equations is made easier if 
the differential equations are expressed in dimensionless form. This 
scales the independent and dependent variables and prevents insignificant 
error terms from dominating the solution. Rewriting equations (1), (2), 

(3), (8), (9), and (10) in dimensionless form yields 


Mm 


d0 

dt^ 


N /a \ 

V ■* 

®\^w/ ' 

e=l 


sin (a + e) - m*W* cos 0 + (11) 

e e o 2 ^ 


★ 

m 


dM 

dt 



(a + Eg) 



- ni*W* sin 0 - TFgH(-z ) 


( 12 ) 


♦ it * 1 2 

= in - 

g o 2 L 



dm 

dt* 



• * 

mf 


(13) 


(14) 
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( 15 ) 


• * • * . ★ 

m = m + 
e a f 


where 


dx 

dt"* 


= M cos 0 


dz 

dt"* 


= -M sin 9 


= 


p c^A 
a w 


t = 




n>^ 


* Pa^. 

X = X 


Pa'S, 


P ■ ~"'o - W 

IIIq 


(16) 

(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


* 

m 


m/m^ 


( 22 ) 


and H is the Heaviside function 


H(s) 



(s ^ 0) 
(s < 0) 


(23) 


Other symbols in equations (11) through (17) have been previously defined. 
Now, all required outputs for the flight module can be computed. The air- 
craft body axis Euler angle 0j^ is related to the angle of attack and 
flight-path incidence angle by 


0j^ = a + 0 

and the wind axis Euler angle relative to the body axis is 


(24) 


0 = -a 


(25) 



The Mach number^ power setting, flap setting, and landing-gear position 
must be expressed in terms of the source time set, which is the input 
control time set augmented with each time when the velocity changes 
by AV. The quantity Av is a user-supplied parameter. In addition, the 
atmospheric properties are expressed as a function of the source time set. 


REFERENCES 

1. Etkin, Bernard: Dynamics of Atmospheric Flight, John Wiley & Sons, 

Inc., C.1972. 

2. Dommasch, Daniel O.; Sherby, Sydney S.; and Connolly, Thomas F.: 

Airplane Aerodynamics, Fourth ed. Pitman Pub. Corp., 1967. 
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TABLE !•- DEFAULT VALUES OF INPUT PARAMETERS 


Input 

parameter 

Default 

Av, m/s 

30 

^ 2 

^w' 

100 

Ag, m2 

IT/4 

h, TC[ 

20 

n>o» kg 

10 000 

N 

1 

ho' " 

0.0 

T 

0.01 

Ee' 

0.0 


10 000 

s 

0.0 

V^, m/s 

0.0 

Xj, m 

0.0 

Yl/ m 

0.0 

m 

0.0 

6l, deg 

0.0 

tn' s 

100 

Vjj, m/s 

125 

m 

10 000 

hn' "> 

10 000 
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Figure 2.- Forces acting on aircraft during flight. 




Lift Coefficient 



Figure 3.- Typical lift coefficient for low-speed 
operation with no ground effect. 
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Drag Coefficient, 



Figure 4.- Typical drag coefficient for low-speed 
operation with no ground effect. 
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Lift Coefficient, 


.8 



Figure 5.- Typical lift coefficient for ground effect. 
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Drag Coefficient 
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.010 
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Landing-gear drag coefficient. 



JTO 


2.4 Jet Takeoff (JTO) Module 
John Rawls. Jr. 

Lockheed Engineering & Sciences Company 


Introduction 

The purpose of the Jet Takeoff (JTO) Module is to calculate the position of an aircraft 
during takeoff. The basic takeoff profile consists of ground roll and climb as shown in 
figure 1. Two optional maneuvers may be api>ende<l to the basic takrsoff proHle. One is 
cutback (in engine power) which is a procedure used to reduce noise levels on the ground. 
The other is a steady turn which may be initiated after a steady-state solution to the 
equations of motion has been obtained. 

An arbitrary fight profile requires a solution to nine differential equations: three force 
equations, three moment equations, and three equations to determine the position of the 
aircraft relative to an Earth-fixed coordinate .system. During takeoff, these equations 
conveniently reduce to four first-order nonlinear differential equations which are solved 
numerically with a fourth-order Runge-Kutta technique. Several assumptions are made in 
this module to simplify the analysis. These assumptions include zero wind, a level runway, 
zero aerodynamic ground effects, and zero weight reduction from the burning of fuel. 

Symbols 


v4||; 

Cd 

<^O.Lg 

Cl 

c 

D 

Fy 

G 

9r 

H 

ha 

L 

^Lg 

M 

Moc 

m 

^eng 

Pw 


aircraft wing reference area, in^ (ft^) 
aerodynamic drag coefficient. | — Bo~r~ 

2 Pn ^ ^ tv 

aerodynamic drag coefficient due to landing gear 
aerodynamic lift coefficient, t — ttt" 

- Air 

speed of sound, m/s (ft/s) 
aerodynamic drag, N (lb) 
ground force, N (lb) 
steady-state climb function 

gravitational constant, 9.8066 rn/s^ (32.1741 ft/s^) 
altitude, m (ft) 

absolute humidity, percent mole fraction 
aerodynamic lift, N (lb) 
landing gear position, Up or Down 
molecular weight of dry air, 28.9644 
aircraft Mach number, V/c 
aircraft mass, W/gr, kg (slugs) 
number of engines 
roll rate, deg/s 
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9w 

R 


r-w 

Te 

Tr 

t 

At 

V 

W 

X 

y 

z 

a 

0 

e 

^tol 

c 

n 

p 

T 

(V'6. ®6> M 

(V^IV, Oyj, 4>w) 

(V’ti»6* ^wbi ^wb) 
n 


pitch rate, deg/s 

univcnial gas constant, 8314.32 iii'“^/K-s^ (49718.96 ft^/°R-s^) 
yaw rate, deg/s 
thrust, N (lb) 

standard sea level temperature. 288.15 K (518.67°R) 
time, s 

source time, s 
incremental time step, s 
aircraft velocity, m/s (ft/s) 
aircraft weight, N (lb) 

mrcraft longitudinal distance from origin, m (ft) 

aircraft lateral distance froni origin, m (ft) 

aircraft altitude above runway, m (ft) 

aircraft angle of attack, deg 

sideslip angle, deg 

Hap control variable, deg 

engine inclination angle, deg 

error tolerance 

dynamic viscosity, kg/m-s (slugs/ft-s) 
integration step size, s 

engine power setting, percent maximum thrust 

air density, kg/m^ (slugs/ft^) 

coefficient of rolling friction 

body axis Euler angles, deg 

wind axis Euler angles, deg 

Euler angles relative to body axis, deg 

turn rate, deg/s 


Subscripts: 

a 

b 

bank 

climb 

cutback 

t 

max 


ambient 
body axis 
steady turn 

climb segment of takeoff 
cutback segment of takeoff 
initial 
maximum 
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min 

rot 

w 

Superscripts: 

* 


minimum 
rotation 
wind axis 

nondimensiunal 

derivative with respect to time 


Input 




The JTO Module requires an extensive set of parameters and tables to define the 
performance characteristics of the aircraft and the engines, the aircraft geometry, and the 
takeoff procedure. The Aircraft Configuration Parameters describe the geometric properties 
of the aircraft. Note in this set of parameters that the wf'ight of the aircraft is the maximum 
takeoff weight since weight reduction due to the burning of fuel is neglected. Aircraft 
and engine performance characteristics are input through the Aerodynamic Lift and Drag 
Coefficient Table, the Landing Gear Drag Coefficient Table, and the Engine Performance 
Table. 

The takeoff procedure begins at brake release and ends when the aircraft reaches a 
designated altitude, a designated distance from brake release, or a designated elapsed flight 
time. The basic takeoff procedure is defined by setting the rotation speed, the desired climb 
speed, and the desired climb angle. In order for the aircraft to become airborne, the speed 
at rotation must be greater than the stall speed and the climb velocity should be at least 
20 percent greater than the stall speed. Note that the climb speed and the climb angle 
are labeled “desired.” If one or both of these parameters is chosen too large, there will be 
insufficient thrust to obtain the anticipated flight profile. Should this occur, reduce 9^ to 
the minimum acceptable climb angle or set 9w to the default value of 2.3*^. This establishes 
an initial climb angle 9w which allows the aircraft to attain the climb velocity as soon 
after liftoff as possible. Once the climb speed is attained, the JTO Module automatically 
computes a new climb angle which achieves the maximum rate of climb for a steady-state 
solution. The maximum angle of attack during rotation ensures that the tail of the aircraft 
does not scrape the runway. 

The takeoff procedure may include two optional maneuvers: cutback and steady turn. 
The cutback procedure is implemented by setting the CUTBACK flag to TRUE, designating 
the altitude at which cutback is to occur, indicating the cutback climb angle, and indicating 
the time required for the cutback procedure to be completed. To include a turn in the takeoff 
procedure, set the BANK parameter to TRUE, indicate the altitude at which the turn is to 
begin, indicate the turn rate, and specify the new flight-path heading. Note that a steady 
turn will not be executed unless a steady-state condition has been achieved. 

The final steps are to establish the origin of the Earth-fixed coordinate system, position 
the aircraft on the runway, and define the initial conditions. Usually, the origin of the 
Earth-fixed coordinate system coincides with the location of the aircraft at brake release; 
therefore, X{ is zero. The initial altitude Z{ indicates the height of the aircraft above the 
runway and should be greater than zero. Setting zi to zero may cause an error in later 
noise calculations if the location of the aircraft and an observer coincide. (This results is 
a singularity in the noise calculation.) The user must also configure the aircraft for takeoff 
by defining an initial angle of attack and an initial flap setting. 

The last input data are the Differential Equation Parameters and the Atmospheric 
Properties Table. Under most circumstances, the default values for the Differential 
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Equation Parameters are adequate. Atmospheric parameters required for the calculation 
of lift, drag, and sound speed are described by the Atmospheric Properties Table. 

Aircraft Configuration Parameters 

aircraft wing reference area, (ft^) 
number of engines 
aircraft weight, N (lb) 
engine inclination angle, deg 

Initial Condition Parameters 

initial time, s 

Xj initial longitudinal distance from origin of Earth- fixed coordinate system, 

m (ft) 

initial altitude above runway, rn (ft) 
ai initial angle of attack, deg 

6i initial flap setting, deg 

Final Condition Parameters 
maximum time, s 

maximum distance from origin, m (ft) 
maximum altitude, m (ft) 

Aircraft Performance Parameters 

desired climb velocity, m/s (ft/s) 
rotation velocity, m/s (ft/s) 
desired flight-path angle during climb, deg 
engine power setting, percent maximum net thrust 
rotation rate, deg/s 

maximum angle of attack during rotation, deg 
coefficient of rolling friction 

Steady Turn Parameters 


BANK 

turn flag 

^bank 

turn altitude, m (ft) 

il 

turn rate, deg/s 

VWk 

desired flight-path heading after completion of turn, Vl>ank positive 
initiates right turn and tlfhaak negative initiates left turn, deg 
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^climb 

n 

a 

Q!max 

r 


^max 

^max 

^niax 


Axij 

^eng 

w 
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CUTBACK 

^cutback 

^cutback 

^^cutback 


At 

^tol 

^max 

^riiin 

a 


Cl 

CD,Lg{Ci) 


H* 

c*{H*) 

p*{H*) 

ha{Hn 


Cutback Parameters 

M 

cutbad^ag “ 

cutbad^altitude, m (ft) 

m 

cutbaclOt'linib angle, deg 
time required to complete cutback, s 
■? 

^ Differential Equation Paraineters 

flight tSfiie increment, s 

error iterance 

*» 

integration step size, s 
maxinattn integration step size, s 
minimi^ integration step size, s 

iJerodynamic Lift and Drag Coefficient Table 
*n 

a 

aircraft^ngle of attack, deg 
flap coatrol variable, deg 

'• r 

aerodystamic lift coefficient, t — 77T~r~ 

I ^PaV^Atr 

aerodvJamic drag coefficient, t — ^ rr“ 

*' 4 iPaV^Atr 

Landing Gear Drag Coefficient Table 


aerodytamic lift coefficient, -j — ■j .j- 

aerodynamic drag coefficient due to landing gear extension 


Atmospheric Properties Table (ATM) 
‘rt' 

altitude, re 

speed (|f sound, re Ca 

air density, re pa 

dynamic viscosity, re pa 

absolute humidity, percent mole fraction 


Engine Performance Table 


n engine power setting, percent maximum net thrust 

Mcx> aircraft idach number 

Tc(n, Moo) net thr'ist per engine, N (lb) 
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Output 

Two output tables are created by the JTO module. The Flight-Path Table gives the 
aircraft ground coordinates and Euler angles in both the body axis and the wind axis 
coordinates systems. The Source Variables Table is a function of source time with eight 
dependent variables including Mach number, engine power setting, flap setting, landing 
gear position, and ambient atmospheric conditions. A new source time is added to the 
Source Variables Table whenever one of the eight dependent quantities changes value, 

Flight-Path Table 


flight time, s 

aircraft ground coordinates, m (ft) 
aircraft body axis Euler angles, deg 

aircraft wind axis Euler angles relative to the body axis, deg 
Source Variables Table 


ts 

source time, s 


Mach number 

n{ts) 

engine power setting, percent maximum thrust 

^/(*») 

flap setting, deg 


landing gear position, Up or Down 

Pai^a) 

ambient air density, kg/m^ (slugs/ft^) 

Ca(t») 

ambient speed of sound, rn/s (ft/s) 

Paita) 

ambient dynamic viscosity, kg/rii-s (slugs/ft-s) 


absolute humidity, percent mole fraction 


t 

{0,9rvb(t),0) 


Method 

The JTO Module defines the takeoff profile of an aircraft relative to a fixed position on 
Earth. Three frames of reference are used to describe the motion of the aircraft. One is 
fiveri to the Earth with the origin placed at brake release, as shown in figure 1. The x axis 
is parallel to the runway with positive x in the direction of takeoff. The y and z axes form 
a right-handed coordinate system with the z axis pointing positive “downward.” The other 
two reference frames are fixed to the aircraft with the origin located at the center of mass. 

The two aircraft reference frames used in this module are shown in figure 2. One reference 
fr ame is the body axis coordinate system denoted by (x{,,y 5 ,Zft). The positive Xf, axis 
extends forward from the center of the aircraft. The yf, and zj, axes form a right-handed 
coordinate system with the zj, axis pointing positive downward. The other reference frame 
is the wind axis coordinate system denoted by (xw,yw,zw)- In the wind axis coordinate 
system, the x«, axis is aligned with the aircraft velocity vector. The wind axis coordinate 
system is used to solve the equations of motion. 

The orientation of the Earth-fixed axes emd the aircraft body axes are defined by the 
Euler an gles denoted by (V’6, 4>b) shown in figure 3. In the wind axis coordinate system, 
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the Euler angles are denoted by itpw,0w,<t>w)^ The body and wind axis 
related by the angle of attack a and the sideslip angle 0 as follows: 

Euler angles are 

V’fe = + 0 

(1) 


(2) 

II 

(3) 

In this analysis, the sideslip angle ^ is always zero. 

The JTO Module computes the Earth-fi.xed coordinates (x,y,z), the body axis Euler 
angles (i>b, 0b, f(>b} and the wind ^lxis Euler angles relative to the body axes (V’u»6>^u;6> <^wb)- 
The wind axis Euler angles relative to the body axes are given by 

1 

il 

1 

II 

(4) 

^wb = 'l>w-'i’b = 0 

(5) 

o 

II 

1 

B 

II 

B 

(6) 


The equations governing the position and velocity of the aircraft in the wind axis coordinate 
system in a still atmosphere (refs. 1 and 2) are 


dV ^ 

m— = > Tf cos (a + e) — — D — sin^u, 

dt ' 
e=l 


(7) 


mVru: = W cos sin <t>w (8) 

mVqtv — 53 sin (a + £) -h Fy *f L - IV cos cos <pyj (9) 

e=l 


I /j 

^ ^ 
at 

Qw = — ^ COS <t>w + Sin dw 

at 

d\j)w ^ , dOyy 

Tw — — cos dyj COS ((>w -JT 

at dt 

dx 

— = V COS9yjCO^\l)w 
dt 

dy 

— = V cosOwsinxI^w 

dt 

dz 

— = --Vsine^ 
dt 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 


Equations (7), (8), and (9) represent the equations of motion of the aircraft in the x^, 
yw, and Zw directions. Equations (10), (11), and (12) are expressions for the roll, pitch, and 
yaw rates, respectively. Equations (13), (14), and (15) give the change in the ground coordi- 
nates as a function of time. These equations are simplified during takeoff since the bsink 
and flight-path heading angles are zero (i.e., <f>^ = 0 and il>w = 0). With this simplification, 
the roll and yaw rates are also zero (i.e., pw = 0 and = 0) and the pitch rate reduces to 


= 


IT 


(16) 
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The system of equations now reduces to the following first-order nonlinear differential 
equations 


dt 

d9w 

dt 


'N. 


eng 


Te cos (a -I- £■) - tF,j - D - W smGu 

e=l 




euR 


^ " Tf^ sin {a s) F(f L ^ W vosOir 
dx 


e=l 


m 

1 


dt 


= V cos Oir 


— = -l/sin^„. 
dt 


(17) 

(18) 

(19) 

( 20 ) 


Equations (17) through (20) are solved numerically witli a fourth-order Runge-Kutta 
technique. A solution is obtained at each flight time specified l>y 

f|iew ~ ^old 'b (21) 

An integration step size ^ is chosen such that 

^ < At (22) 

and 

^ < ^max (23) 

The integration scheme adjusts the step size to meet the desired error tolerance. Figures 4 
and 5 show the forces acting on the aircraft during the ground roll and climb. The ground 
force term Fg in equations (17) and (18) represents the resistance due to friction resulting 
from contact between the aircraft wheels and the runway. This force is positive as long 
as the wheels remain in contact with the runway and Ix'comes zero at the point of liftoff. 
An expression for the ground force is obtained by noting that the flight-path angle and 
pitch rate K in equation (18) are zero during ground roll. Consequently, 





'^fllg 

L - ^ TeSin(a -I- -:) 
e=l 


[z =0) 

( 2 < 0 ) 


(24) 


The rolling friction coefficient r in equation (17) is a function of the landing gear and the 
surface characteristics of the runway. Two assumptions are made to simplify this parameter. 
The surface is assumed to be uniform during the ground roll and the ground force Fg is 
assumed, to result solely from the main landing gear. This assumption allows the friction 
coefficient to remain constant until liftoff. The aircraft lift and drag are computed by 


L=\paV‘^Ar,[Ci{cc.bf)] 


(25) 


and 


D=l 


\po (Oi^/) + <7o,Lg (Cl,)] 




\paV'^A^ [Cd (a, 6;) + cos 

Cp,Lg(Ct)| 

(3 S ^^climb — 

(26) 

\paV^A^CD (a, 6/) 
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In equation (26), At riimh is 

^^cUmb — ^ ~ ^clinib (27) 

and tciimb is the time at the start of the climb stage. The lift coefficient Cl and the drag 
coefficient Cp are functions of the aircraft angle of attack and flap settings. An additional 
source of drag Cd , which is a function of the lift coefficient, is due to the extension of 
the landing gear. Landing gear drag is present during the ground roll stage and the first 
3 seconds after liftoff. The landing gear drag coefficient is multiplied by a cosine term 
to ensure a smooth transition in the drag force. These coefficients are obtained from the 
Aerodynamic Lift and Drag Coefficient Table and the Landing Gear Drag Coefficient Table. 

The net thrust Te is a function of the engine power setting and the flight Mach number. 
Values for the net thrust are obtained from the Engine Performance Table. 

Takeoff Procedure 

The takeoff procedure is divided into four stages: ground roll, climb, cutback, and steady 
turn. Ground roll and climb are the basis for all takeoff profiles. Cutback and steady turn 
are optional stages. 

Ground ivll. Ground roll begins with the aircraft at rest. During the ground roll stage, 
the climb angle 9u.' and the pitch rate &tv are zero. With these restrictions, the problem is 
one-dimensional and requires solutions to equations (17) and (19). The initial conditions 
for the ground roll stage arc 

V(ti) = 0 (28) 

x{ti) = xi (29) 

The angle of attack a, climb angle 6w, and coordinates y and z remain constant during 
ground roll; that is, a = q^, = 0, y = 0, and z = z,. Equations (17) and (19) are solved 

iteratively until the rotation velocity is achieved. After reaching the rotation velocity, the 
angle of attack is increased by 

anew = oiold + « At (30) 

until the maximum angle of attack has been achieved. As the nose of the aircraft rotates 
“upward,” the velocity and the angle of attack continue to increase. The velocity vector 
however remains parallel to the runway. The ground roll stage terminates when the ground 
force term Fg becomes zero. 

Climb. Climb begins at the point of liftoff when the aircraft main landing gear 
leaves the ground, and 3 seconds into the climb stage, the landing gear is automatically 
retracted. During the climb stage, the aircraft accelerates to the desired climb velocity, 
and simultaneously, the climb angle increases from zero until the desired climb angle is 
obtained. The aircraft coordinates are obtain^ by solving equations (17) through (20). 
During the solution of these equations, one of the following conditions will arise: 

Condition 1 ( V < Vdimb and 6^, < ^dimb): Eyery climb stage begins with the aircraft 
velocity and climb angle less than the desired values. The angle of attack is evaluated 
at the beginning of the climb stage to determine if rotation has been completed. If 
rotation has not been completed, the angle of attack continues to increase until the 
maximum angle of attack is achieved or until the aircraft acceleration rate becomes 
negative. If the acceleration rate becomes negative, the aircraft has “overrotated” ; 
in which case, the angle of attack is reduced by 1 deg/s until the acceleration rate 
becomes positive. Equations (17) through (20) are solved numerically until one of the 
conditions described subsequently occurs or until a steady-state solution is obtained. 
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Condition 2 (V = Kdimb < ^dimb): aircraft attains the climb 

velocity before attaining the desired climb angle, the angle of attack is adjusted so 
that the aircraft no longer accelerates. The appropriate angle of attack is computed 
by setting equation (17) to zero 

^eriR 

T(,cos(q + £) - D - U sinO, 

e=l 

and solving for a numerically. The remaining energy in the system is used to gain 
altitude. If additional energy is available after tlu' climb angle is attained, a second 
climb segment is automatically initiated as descrilx'd by condition 4. 

Condition 3 (K < K-limb an<l 0,^ = ^dimb): The alternative to condition 2 is when 
the aircraft attains the climb angle before attaining the desired climb velocity. For 
most takeoff procedures, setting a climb angle less than the maximum climb angle 
attainable by the aircraft is desirable. This establishes ^,.u,nb as an initial climb 
angle which is held constant to allow the climb velocity to be achieved as quickly as 
possible. To ensure that the initial climb angle is maintained, equation (18) is .set to 
zero and solved numerically for q: 


= 0 


(31) 


y"* Tf. sin {a + e) + L 

C=1 


ir cos 6 , 1 , 


= 0 


(32) 


The angle of attack now controls the climb angle, and the remaining energy in 
the system is used to accelerate the aircraft. If additional energy is available after 
the climb velocity is attained, a second climb segment is automatically initiated as 
described by condition 4. 

Condition 4 (K = and 9w = ^dimb): Condition 4 arises when the aircraft 

is capable of a greater climb speed or climb angle than is specified by the input 
parameters, ^^hen this condition occurs. Iciimb held constant and the excess 
energy is used to increase the rate of climb. .A new climb angle is computed by 
setting the right side of equations (17) and ( 18) to zero and simultaneously solving 
for 6 w The new climb angle can be written in terms of the lift-drag ratio and the 
thrust-weight ratio as 


Oy, = sin 


f {L/D)G 


■ {L/D)G ■ 

2 

1-G2 

- > 

■ 

\ 1 + {L/Df \ 


1 + {L/D)\ 


1 + {L/Df 



(33) 


where the steady-state climb function G is defined as 


N, 


eng 


G = 


E 

e=l 


Te . 


w 


— cos (tt -I- =) 4- sin (a -t- e) 


(34) 


Steady-state solution: A steady-state solution is defined when both the acceleration 
rate and the pitch rate satisfy the following criteria: 


m 


dt 


^ ‘tol 


(35) 
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and 


mV 


dt 


< £tol 


The Earth-fixed coordinates can now be computed by 


^new ^old "i* ^ ^tcos6^jt* 


(36) 


(37) 


and 

Zncw = 2 old “ ^ (^®) 
The climb stage ends when one of the following conditions are met: 

X ^ XiiiHx (39) 

Z ^ 'Ziiinx (^9) 


or 


t ^ ^niax 


(41) 


The climb stage may also end when the cutbfick altitude is obtained or when a 
steady-state solution is obtained and a steady turn maneuver is to be performed. 


Cutback. Cutback is initiated when the aircraft reaches the user-specified cutback 
altitude. With the aircraft angle of attack and the climb velocity remaining constant 
throughout the maneuver, the thrust required to obtain the cutback climb angle is computed 


i>y 


.^cutback 


W [cos ^cutback i^/^) ^cutback] 
(L/D) cos (a -(- £) -t- sin (« -I- e) 


(42) 


Thrust is reduced linearly such that 


^eiig 

-Veiig 


^ ^ ^C,IICW 

= ^ 




C = 1 





where the thrust increment given by 

e=l 


^enR 

£ AT, 


C=1 


(^cng 


Tft (n) ^cutback I 
^^cutback 


) 


(43) 


(44) 


and rfi(n) is a function of the power setting at the beginning of the cutback maneuver. The 
climb angle corresponding to each new thrust setting is computed by equation (33) with 
Te = Te,new in equation (34). At the end of Afcutback. in seconds, cutback is completed 
and a steady-state solution to equations (17) and (18) is obtained with the new climb angle 
equal to the cutback climb angle. 

Steady turn. When an aircraft executes a steady turn, the aircraft velocity, climb 
angle, and roll angle are required to be constant. Consequently, a steady-state solution 
to equations (17) and (18) is required before a turn can be initiated. Once this condition 
is met, a steady turn is initiated as soon as the user-specified turn altitude is achieved. 
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The Hight-path heading during the turn varies as 

— V'uJ,old ^ ('^^) 

where the turn rate Q is constant. A positive value for fl produces a right turn, from the 
pilot’s perspective. Under these oanditions, the yaw rate is 

>„, = fjcos0,t COS</<u; (46) 

and the balance of forces in the direction is 

mVQaas 6,^ cos <p„, — W cos 0^. sin <j>„- (47) 


when solving equation (47) for an angle 

-1 //1a^ 

j <p^, = tan ^ (48) 

Hr 

was found to be required to exeoite the turn. The Earth-fixed coordinates are obtained 
from the following equations: 

x„ew = x„M -H cos 8^. [sin ~ sin (t/;,, - il At)] (49) 

ynew = yold + cos 0w [cos (ti'u, - n AT) - cos ri^w] (50) 

^new — ^old “ sin 6*]^ (51) 

The steady turn procedure is terminated when 

Ww — l/^baiik (52) 

or when x, z, or t exceed the maximum values given by equations (39), (40), and (41). 
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2.5 Jet. Landing (JLD) Module 
MarlcWilson 

Lockheed Engineering & Sciences Company 

Introduction 

The Jet Landing (JLD) Module computes the position of an aircraft during an approach 
to the runway based on the aircraft performance characteristics. The landing procedure can 
be divided into as many as five segments. The basic landing profile consists of one segment 
with an approach angle of 3°. During each segment, the approach speed and flight-path 
angle are required to be constant. This constraint reduces the differential equations of 
motion to an algebraic form which is used to determine the angle of attack and thrust 
required to maintain the input approach speed and flight-path angle. 

Symbols 


Cd 


Cl 

c 

D 

Fo 

9r 

H 

ha 

L 

M 

Moo 

m 

^eng 

Naeg 

R 

Te 

Tr 

t 

ia 


aircraft wing reference area, m"^ (ft^) 
aerodynamic drag coefficient. -> — — 

^p„V-Au’ 

aerodynamic drag coefficient due to landing gear 

aerodynamic lift coefRcient, -i — ^ — 

^Pay^Au- 

speed of sound, m/s (ft/s) 
aerodynamic drag, N (lb) 
ground force, N (lb) 

gravitational constant, 9.8066 in/sP (32.1741 ft/s^) 
altitude, m (ft) 

absolute humidity, percent mole fraction 

aerodynamic lift, N (lb) 

landing gear position. Up or Down 

molecular weight of dry air, 28.9644 

aircraft Mach number, V/c 

aircraft mass, Wlg^., kg (slugs) 

number of engines 

number of flight-path s^ments 

universal gas constant, 8314.32 m^/K-s^ (49718.96 ft^/°R-s^) 
thrust, N (lb) 

standard sea level temperature, 288.15 K (518.67**R) 
time, s 

source time, s 
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iucrc-rnental time step, s 
aircraft velocity, m/s (ft/s) 
aircraft weight, N (lb) 

aircraft longitudinal distance from origin, m (ft) 

aircraift lateral distance from origin, m (ft) 

aircraft altitude above runway, m (ft) 

aurcratft angle of attawik, deg 

flap control variable, deg 

engine inclination amgle, deg 

dynamic viscosity, kg/m-s (slugs/ft-s) 

engine power setting, percent maximum thrust 

thrust reverse power setting, percent maximum thrust 

air density, kg/m^ (slugs/ ft^) 

coefficient of rolling friction 



body axis Euler angles, deg 


wind axis Euler angles, deg 

(V^w6> 

Euler angles relative to body axis, deg 

Subscripts: 


a 

ambient 

b 

body axis 

f 

final 

H 

runway threshold 

i 

initial 

Lg 

landing gear 

n 

counter 

w 

wind axis 

1 

first segment 

2 

second segment 

Superscript: 


* 

nondimensional 

Input 



A description of the aircraft geometry is required. The Aerodynamic Lift and Drag 
Coefficient Table and an Engine Performance Table describe the characteristics of the 
airframe and engine. The atmosphere is described by the Atmospheric Properties Table. 
Input parameters are used to define the number of landing segments and other altitude 
dependent variables such as the altitude for landing gear extension. An array of altitudes. 


At 

V 

W 

X 

y 

Z 

a 

h 

e 

n 

nRv 

p 
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flight-path angles, and approach velocities for each segment of the landing profile 
in the Landing Profile Table. 

Input Parameters 


Nfieg 

number of segments 

IIrv 

thrust reverser power setting, percent maximum net thrust 

At 

time increment, s 


final aircraft velocity, m/s (ft/s) 

W 

aircraft weight, N (lb) 

^Lg 

altitude for flap and landing gear extension, m (ft) 


altitude for end of runway crossing, m (ft) 

h 

landing configuration flap setting, deg 


Landing Profile Table 


array of approach altitudes, m (ft) 

^w,i 

array of approach flight-path angles, deg 
array of approach velocities, m/s (ft/s) 


Aircraft Configuration Parameters 

Au’ 

aircraft wing reference area, m^ (ft^) 

■^cng 

number of engines 

£• 

engine inclination angle, deg 

T 

coefficient of rolling friction 


Aerodynamic Lift and Drag Coefficient Table 

a 

aircraft amgle of attack, deg 

h 

flap control variable, deg 

CL(a,6f) 

aerodynamic lift coefficient, 


aerodynamic drag coefficient, -j — , a 


Laoiding Gear Drag Coefficient Table 

Cl 

aerodynamic lift coefficient, ^2^," 

Cdm^Cl) 

aerodynamic drag coefficient due to landing gear extension 


are defined 
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Atmospheric Properties Taijlc ( ATM) 


H* 

altitud^, rej^ 


speed aif sound, re Ca 

p*(H*) 

air desity, re pa 


dynamic viscosity, re pa 

ha{H*) 

absolute humidity, percent mole fraction 


Engine Performance Table 

n 

engine-power setting, percent maxiimim net thrust 

Moo 

aircraft Mach number 

Te{n,Moo) 

net thmst per engine, N (lb) 


Output 

Two output tables are created by this module. The Flight-Path Table gives the aircraft 
ground coordinates aiul Euler angles in both the body axis and the wind axis coordinates 
systems. The Source Variables Table is a function of source time with eight dependent 
variables including Mach number, engine power setting, flap setting, landing gear position, 
and ambient atmospheric conditions. A new source time is added to the Source Variables 
Table whenever one oFthe eight dependent quantities changes value. 

> Flight-Path Table 


t 

(l(t), 0,2(0) 
( 0 ,<> 6 ( 0 , 0 ) 


A/oo(0) 

n(o) 

6f{ts) 

^Lg(*a) 

Pa{ta) 

Co(0) 

ha{ts) 


flight time, s 

■; aircraft ground coordinates, m (ft) 

* aircraft body axis Euler angles, deg 
aircraft wind axis Euler angles relative to body axis, deg 

' Source Variables Table 

I source time, s 
Mach number 

engine power setting, percent maximum thrust 
flap setting, deg 

landing gear position. Up or Down 
ambient air density, kg/m^ (slugs/ft^) 
ambient speed of sound, m/s (ft/s) 
ambient dynamic viscosity, kg/m-s (slugs/ft-s) 
absolute humidity, percent mole fraction 
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Method 

The flight trajectory generated by this module is based on balancing the forces in the 
horizontal direction by ~ 


IVeng jy 

-Tpg + ^2 COS (a + e) - \V sin^u> - D = (1) 

C=1 

and in the vertical direction by 

^'eng ^ 

Fj + ^ Tfl sin (o + £) - cos0,y + L — mV (2) 

C=1 

These equations (refs. 1 and 2) are defined in the wind axis coordinate system shown in 
figure 1. The ground force term Fg is defined by 


{ 'Veng 

W - L - 53 

e=l 

0 

and remains zero until the moment of touchdown. The coefficient of rolling friction t is a 
function of landing gear and surface characteristics. 

The aircraft lift L and drag D are computed by 

l = \paV^A,„[CL{a,6f)] (4) 


(2 = 0 ) 
( 2 > 0 ) 


(3) 


and 

D = ^paV^A^ [Cd (a, 8f) + Cd, L g [Cl)] (5) 

where the lift and drag coefficients are obtained from the Aerodynamic Lift and Drag 
Coefficient Table and the Landing Gear Drag Coefficient Table. 

The position of the aircraft as a function of time is given by 

^ = Vcos6w ( 6 ) 

at 

and , 

^ = V'sinO,. (7) 

This module divides the landing procedure into segments defined by the parameter iVseg- 
An array of altitudes, flight-path angles, and approach velocities are specified in the Landing 
Profile Table. The arrays are defined such that the first value corresponds to the outermost 
segment. Figure 2 shows an example of a two-segment landing profile. Note that the flight- 
F>ath angles are defined to be negative down. The origin of x is defined at the end of the 
runway. The initial distance of the aircraft from the end of the runway is computed from 
the array of altitudes and flight-path angles by 





n=l 


t&n 


( 8 ) 
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where the last altitude in the summation corresponds to the altitude at the end of the 
runway; that is, 

•*iVseg+l = 

By maintaining a constant approach speed and flight-path angle, the derivatives in 
equations (1) and (2) are zero. This allows the use of the following algebraic equations 
(assuming the aircraft is not on the ground; i.e., Fg = 0); 

^eng 

^T,cos{a + e)-Wsmd„.- D = 0 (10) 

e=l 


and 


Te sin (q -h e) — Vr cos + L = 0 (11) 

e=l 


To calculate the aircraft angle of attack a required to maintain the desired approach speed, 
a small angle approximation is made for a -I- £ in equation (11). Assume that sin (q -I- c) 0 
and equation (11) becomes 

L = Wcos6u, (12) 

which can be written in terms of the lift coefficient as 


_ W cos $w 
^ ~ ^aV^Au, 


(13) 


The angle of attack that satisfies equation (13) is found by interpolating the Aerodynamic 
Lift and Drag Coefficient Table. 

The thrust required to maintain the desired flight-path angle is obtained from 




e=l 


W [cos 0w -{-{L/D) sin 0^] 
{LjD) cos (a -t- e) -(- sin (a + e) 


(14) 


During the airborne segments of the landing trajectory, the approach velocity remains 
constant so that the aircraft coordinates are 


Xj = Sj_i -I- Kcos^u; At ( 15 ) 

and 

Zj = Zj-\ - V sin^w At (16) 

where j indicates segment number. After touchdown, the aircraft velocity becomes a 
function of time due to the ground force term. The velocity then becomes 

Vj = Vj-i+AV (17) 


where the reduction in velocity is given by 



m 


-rFg - D - mgr sin Ow+^^'I'e cos (a -t- e) 

e=l 


At 


(18) 
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Thrust Tt is calculated by using the Engine Performance Table where engine power setting 
and aircraft Mach number are used as inputs. After touchdown, thrust reversing may be 
employed to rapidly decelerate the aircraft. The engine power sei is replaced by the 
thrust reverser parameter as input to the Engins Performance Table, i hrust is reversed by 
assigning a negative value to the thrust reverser parameter Orv For example, a value of 
IIrv = -0.5 would provide thrust reversing at a power setting of one half full power. 

The aircraft Euler angles are also required as output (fig. 3). Since the two-degree- 
of-freedom assumption has been made, the landing profile takes place only in the vertical 
plane (i.e., no yaw or roll motion). This assumption implies that the values of the Euler 
angles and (ph are both zero. The aircraft body axis Euler angle Of, can be calculated by 
the following relation between angle of attack ^d flight-path angle: 

Ofj = ot 9x0 ( 19 ) 

The body axis Euler angle 9i is transformed to the wind axis by using the following equation: 

Kh = 


References 

1. Etkin. Bernard: Dynamicn of Atmoxphmc Flight John Wiley Ai Sons. Inc., c.1972. 
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Pitman Publ. Corp., c.1967. ' 
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2.6 Steady Flyover (SFO) Module 

John Rawls, Jr. 


SFO 


Lockheed Engineering & Sciences Company 


Introduction 

The Steady Flyover (SFO) Module defines the position of an aircraft based on kinematics 
(ref. 1). This approach enables flight paths to be defined without requiring engine and 
aircraft performance data as input. Aircraft motion is restricted to constant velocity or 
uniform acceleration along a rectilinear flight path. Flight paths that do not conform to 
this restriction may be divided into segments. A segmented flight path is constructed by 
executing this module once for each segment. 

Input to the SFO Module consists of parameters that define a single flight-path segment. 
Output consists of parameters that allow repeated executions of the module, and two tables 
that are used by other ANOPP (Aircraft Noise Prediction Program^ dules to predict 
noise produced by a moving source. One table defines the position of the aircraft as a 
function of time, and the other identifies changes in the noise source characteristics. 

Symbols 


a 

c 

d 

9r 

H 

h(i 

^Lg 

A/ 

A/oC 

N 

R 

Tr 

t 

^Lg 

At 

V 

X 

y 


aircraft acceleration, m/s^ (ft/s^) 
speed of sound, in/s (ft/s) 
distance, m (ft) 

gravitational constant, 9.8066 in/s^ (32.1741 ft/s*) 
altitude, in (ft) 

absolute humidity, percent mole fraction 
landing gear position. Up or Down 
molecular weight of dry air, 28.9644 
aircraft Mach number, V/c 
number of incremental time steps 

universal gas constant, 8314.32 m*/k-s* (49718.96 ft*/°R-s*) 
standard sea level temperature, 288.15 K (518.67*’R) 
time, s 

landing gear reset time, s 

source time, s 

incremental time step, s 

aircraft velocity, m/s (ft/s) 

aircraft longitudinal distance from origin, m (ft) 

aircraft lateral distance from origin, m (ft) 

aircraft altitude above runway, m (ft) 
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a 

aircraft angle of attack, deg 

tf 

flap control variable, deg — 

j 

dynamic viscosity. kg/m*s (slugs/ ft-s) 

n 

engine power setting, percent maximuiii 

p 

air density, kg/nv^ (slugs/ft^) 


body axis Euler angles, deg 

( t’u*. Oic- Ou’ ) 

wind axis Euler angles, deg 

( Vu i)- ^tcl* 

Euler angles relative to body axis, deg 

Subscripts: 

a 

ambient 

b 

body axis 

f 

final 

i 

initial 

ref 

reference 

ir 

wind axis 

1 

first segment 

2 

second segment 

Superscript: 

>;C 

nondimensional 


Input 

\ Bkhl proWe mav cousiat of aa ma..y aesmc.c. aa Che usee desire Each 
A nigm pionie . Module. To create a multiple segment Hight 

"*r',hrAP?END ZIIZ be se° ,o TREE, This cases .he da.a ccea.ed foe Che 

n°gh':P.thTablea,cd,heS^^^^^ 

;rc=T„X':pp;opei«: ::^e.'‘.oSe append pa„, ue.ee is FALSE, each e.ecuuou 
of the SFO Module creates a new flight profile. 

The Mach ...mebee upda.e cei.eeion paeame.ee A.l/x ideu.ifi« evhen .he flight .yeed has 

coucaius au e„,e>- foe 

of “““ ehe uumbee of source uoise peedic.ous 
T: a'°e"relc?Xti?oL;« ..me by se..,„g .he Maeh uumbee upda.e 
criteeiou parametee. Generaliy. A.i/x - O.Oo is adequate. 

The couepman^al 

"anSe^^ : SomeiimL the inclination angle is a more convenient parameter for constructing 
a flight-path segment. 
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The position of an aircraft is calculated with the initial and final condition parameters. 
The Aircraft Configuration Parameters do not affect the flight profile but are required as 
output by the Source Variables Table. Information from the Source Variables Table is 
used by other ANOPP modules to predict the noise characteristics of the engine and the 
airframe. With the exception of the landing gear, the Aircraft Configuration Parameters 
remain constant throughout a flight segment. The Earth-fixed coordinated system is defined 
such that z is positive down as shown in figure 1. For convenience, the z coordinates 
(corresponding to altitude) are input as positive values and converted to negative values 
within the module. In order to compute the aircraft Mach number, the speed of sound is 
provided by the Atmospheric Properties Table. 

Input Parameters 


APPEND 

Ji 

AA/oc 

Af 


ZOPT 


multiple segment flight-path flag, True or False 

initial step number 

Mach number update criterion 

time step increment, s 

altitude of runway above reference level, m (ft) 

computational option flag, for ZOPT = 1, input zj and disregard and 
for ZOPT = 2, input 9w and disregard zj 

inclination angle of flight vector with respect to horizontal, positive for 
climb and negative for descent, deg 


Initial Condition Parameters 


ti initial time, s 

Vi initial velocity, m/s (ft/s) 

X| initial longitudinal position from origin, m (ft) 

yi initial lateral position from origin, m (ft) 

Zi initial altitude above runway, m (ft) 

Final Condition Parameters 


*/ 

final 

Kr 

final 

XJ 

final 

yj 

final 


final 


time, s 

velocity, m/s (ft/s) 

longitudinal distance from origin, m (ft) 
lateral position from origin, m (ft) 
altitude above runway, m (ft) 
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Aircraft Configuration Parameters 


^Lg,t 

initial landing gear position, Up or Down 

*Lg 

landing gear reset time, s 

n 

engine power setting, percent maximum thrust 


flap setting, deg 

a 

angle of attack, deg 


Atmospheric Properties Table (ATM) 

H* 

altitude, re 

c*(if*) 

speed of sound, re Ca 

p*(/T) 

air density, re pa 


dynamic viscosity, re pa 

ha{H*) 

absolute humidity, percent mole fraction 


Output 

The Final Condition Parameters provide pertinent information necessary for a repeat 
execution of the Steady Flyover Module. When a multiple-segment flight profile is created, 
the final conditions are used as the initial conditions for the next segment. Two output 
tables are created by this module. The Flight-Path Table gives the aircraft ground 
coordinates and Euler angles in both the body a.\is and the wind axis coordinates systems. 
The Source Variables Table is a function of source time with eight dependent \ariables 
including Mach number, engine power setting, flap setting landing gear position, and 
ambient atmospheric conditions. A new source time is added to the Source Variables 
Table whenever one of the eight dependent quantities changes value. 

Final Condition Parameters 

If the SFO Module is executed with APPEND = TRUE, the following parameters 
become the input parameters for the next execution of the SFO Module: 

final step number 

final landing gear position. Up or Down 
final time, s 

actual final velocity, m/s (ft/s) 
actual final longitudinal position from origin, m (ft) 
actual final lateral position from origin, m (ft) 
actual final altitude above runway, m (ft) 


^Lg,/ 

Vf 

Xf 

Vf 
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Flight-Pat^ Table 


t flight time, s 

{x{t)^y{t),z{t)) aircraft position coordinates, in (ft) 

aircraft body axis Euler angles, deg 

{O,0^lj{t),O) aircraft wind axis Euler angles relative to body axis, deg 

Source Variables Table 



source time, s 

^oci^s) 

Mach number 

n(^) 

engine power setting, percent maximum thrus 

6j{ts) 

flap setting, deg 


landing gear position. Up or Down 

Pai^s) 

ambient air density, kg/m'^ (slugs/ft^) 


ambient speed of sound, m/s (ft/s) 

Pai^h) 

ambient dynamic viscosity, kg/m-s (slugs/ft-s) 


ambient absolute humidity, percent mole fraction 


Method 

The SFO Module computes the position of an aircraft based on kinematics. This 
approach is useful when aircraft performance data arc not available or when simple flight 
profiles are required. Figures 2 through 4 illustrate three examples of flight profiles that 
can be constructed by using the SFO Module. 

Figure 2 illustrates a single-segment flight path where the velocity and the altitude 
remain constant. A flight-path segment is defined by the initial and final position of the 
aircraft and by the initial and final velocity. At the beginning of each segment, the angle 
of attack, the flap setting, and the power setting are provided as input. As long as the 
flight path can be defined between two points, the motion of the aircraft may be climbing, 
descending, accelerating, or decelerating. The aircraft motion may also be stationary; 
therefore, a flight profile may be created for a hovering aircraft. 

The flight profiles shown in figures 3 and 4 are constructed by repeated executions of the 
module. The landing profile shown in figure 3 is divided into two segments, representing 
approach and ground roll. The final position and velocity for the approach segment becomes 
the initial conditions for the ground roll segment. 

Figure 4 illustrates a takeoff profile that has been divided into four segments: ground 
roll, first climb segment, second climb segment, and cutback. During ground roll and the 
first climb segment, the aircraft is accelerating. The landing gear may be retracted during 
the first climb segment by specifying the new gear position and the reset time which is 
measured from the beginning of the first climb segment. During the last two segments, 
the aircraft velocity is constant, but the two segments differ in climb angle and in power 
setting. 
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The position coordinates of the aircraft are computc<l by the following equations: 


= Ifc + 1 + 


(' 


Zk+l = ( 


At^j 

1 COS Oj,) cos 

(1) 

At") 

1 

(2) 

^u At") sin0,„ 

(3) 


Vk At + 


Equations (1), (2), and (3) are solved iteratively for k = ./ to N, where N is the integer 
part of 


N = Ji + 


At 


(4) 


In equation (4), J, is the initial step number, t is the total time required to complete the 
flight-path segment, and At is the incremental time step. The aircraft velocity at time 
step k is 

Vk = Vk.i-\-nAt (5) 

where V*._i is the velocity from the previous time step and a is the acceleration. The 
acceleration of the aircraft given by 


a = 


Vj - vl 

2d 


( 6 ) 


and is computed from the initial and final velocities and the distance traveled during the 
segment, which is 


d= \J Ax^ -H Ai/" + Az“^ 
The incremental distances Ax, Ay, and Az are 

(7) 

Ax — Xf - j‘, 

(8) 

> 

II 

1 

(9) 

> 

N 

II 

1 

(10) 


The total time required to complete a segment is 


t = 


'tf-ti 
^ 2d 

[vrrvj 


(Vi = Vf=0) 

(|a| > 0) 


( 11 ) 


The flight-path heading and the climb angle 0^ remain constant during a segment. The 
flight-path heading is given by 

ipw = tan~^ ^ (12) 

An option is provided which aJlows either the final altitude Zf at the climb angle to be 
specified as input. If is specified, the climb angle is comput^ by 


0y, = tan 1 


Ax 


(13) 


rev. 11-93 


2 . 6-6 



SFO 


If i> >pecititHi. A: i> computed })v ~ 

A: = A.rtan^y,, ' 14 1 

The aircraft body axis Euler angles and t ^ an^ ( alciilated by using the following 
relationship between angle of attack. Hight-path angle, and Hight-path heading. See tig. o. » 
The eliml) angle in the body axis is 


= M - « 


and the Highi-j)ath heatiing in the body axis is 

L l, = (. V 


. lo- 


- 10 i 


The body axis E\iler angle is transtunned to the wind axis with tla llowing relation: 


^ K = 17 . 

This ino«itile creates two output tal)les. One is the Fliglir-Patli Table, which defines 
the aircraft position and Euler angles for each time increment A^ The other is the 
Source \ariables Table, which identities times during the Hight when the characteristics 
of the source noise an* altered by dianges in engine power setting, aircraft configuration, 
ambient conditions or aircraft velocit\'. The Source \’ariables Table always contains an 
entry correspondinu to the start of a new flight-path .'-egment. 

The flight Mach numl)er. given by 


c// 


is evaluated from the aircraft velocity and the ambient sound speed. The sound speed ci //“ } 
is obtaineil from the Atmospheric Pro])erties Table as a function of the nondimeiisional 
altitude //" given by 


//" = 


RTr Mgr 


( 19 ) 


The parameter is the height of the runway alcove the reference level (usually sea level) 
specified in the Atmospheric Properties Table. 


Reference 


1. Btvr. FvriiiiicUiti P.; hikI .luiinstuii. E. Ru>>t'll. .Ir.! > for utid Dytiotutt s, 

McGraw-Hill Book Cu.. U)G-. 
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3. PROPAGATION EFFECTS 




3.1 ATMOSPHERIC ABSORPTION MODULE 


INTRODUCTION 

As sound waves propagate from the source to the observer, they are 
attenuated due to atmospheric cibsorption. The Atmospheric Absorption 
Module accounts for the decrease in noise intensity produced by the atmo- 
sphere. This sound attenuation is expressed as an absorption coefficient 
which depends on frequency, temperature, pressure, and humidity. 

The temperature, pressure, and humidity are assumed in ANOPP to be 
functions of altitude only. Therefore, the absorption coefficient is 
expressed as a function of altitude and frequency. It is shown that a 
mean absorption coefficient, which is the average value from ground level 
to a given altitude, can be applied to the sound intensity. 

The atmospheric absorption occurs due to four basic causes. First is 
the classical loss due to thermal and viscous effects. Second is the 
molecular loss due to the rotational relaxation of nitrogen and oxygen 
molecules. Third is the molecular loss due to the vibrational relaxation 
of nitrogen molecules. Fourth is the molecular loss due to the vibra- 
tional relaxation of oxygen molecules. The total absorption coefficient 
is the sum of the coefficients for each of the four effects. 

This module computes a table of the atmospheric-absorption coeffi- 
cient in dimensionless form as a function of frequency and dimensionless 
altitude. This table is then available to be used with the Propagation 
Module to correct sound intensities for the atmospheric-absorption effect. 


SYMBOLS 

c speed of sound, m/s (ft/s) 

f frequency, Hz 

f^j^ relaxation frequency, Hz 

9 ^ acceleration due to gravity, m/s^ (ft/s^) 

H altitude, m (ft) 

Hj^ ground level altitude, m (ft) 

h absolute humidity, percent mole fraction 

M molecular weight of air 

P pressure, Pa (Ib/ft^) 
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<p^> mean-square pressure, Pa^ (Ib^/ft^) 

R universal gas constant, m^/K-s^ (ft^/®R-s^) 

r distance, m (ft) 

T temperature, K (^R) 

X fractional molar concentration 

y dimensionless altitude, Mg^(H - H 2 ^)/rT^ 

a absorption coefficient, nepers/m (nepers/ft) 

0 characteristic vibrational temperature, K (°R) 

\x dimensionless absorption coefficient, ac^/f, nepers 

p reference density, kg/m^ (slugs/ft^) 

Subscripts : 
cl classical 

n nitrogen 

o oxygen 

r standard sea level value 

rot rotational 

s source 

vib vibrational 

Superscript: 

* dimensionless value 


INPUT 

The two basic independent variables for the module are frequency and 
altitude. Center-frequency values for 1/3-octave bands and the altitude 
values in the Atmospheric Module (ATM) define the independent-variable 
inputs. In addition, the module requires several of the atmospheric 
properties produced by the ATM. 

f 1/ 3-octave-band center frequency, Hz 
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Atmospheric-Properties Table 
y dimensionless altitude, Mg^(H - 

h(y) humidity, percent mole fraction 

p*(y) pressure, re p^. 

T*(y) temperature, re 


OUTPUT 

The output is a table of the dimensionless absorption coefficient as 
a function of frequency and altitude. 


Absorption-Coefficient Table 
f frequency, Hz 

y dimensionless altitude, Mg^(H - 

y(f,y) dimensionless absorption coefficient, ac^/f, nepers 


METHOD 

The sound intensity lost due to atmospheric absorption is expressed 
as an atmospheric-absorption coefficient. This coefficient is defined by 
the relation 


<P^(r)> 




( 1 ) 


2 

where <^P is the mean-square acoustic pressure at some distance r 

from the source, <p^ (rg)> is the mean-square acoustic pressure at the 

source, and a is the absorption coefficient. The absorption coefficient 
is expressed as a sum of four components as 


a 


a 


cl 


“rot ^ 


a . , + 

vib,o 


a. 


vib,n 


(2) 


where is the classical loss due to thermal and viscous effects, 

a is the molecular-absorption loss due to the rotational relaxation of 

rot 

oxygen and nitrogen molecules, ot . , is the molecular-absorption loss 
due to the vibrational relaxation of oxygen molecules, and is 

the molecular-absorption loss due to the vibrational relaxation of nitro- 
gen molecules. Additional sources of atmospheric attenuation are assumed 
to be negligible. 
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As a result of extensive theoretical and experimental effort, 
expressions have been developed for each of the terms in equation (2) . 
Sutherland (ref. 1) presents the theoretical development and all existing 
experimental data. Reference 2 further refines that work. All of the 
empirical equations used in this module are results of the work in these 
references . 

An important parameter in the vibrational-relaxation absorption 
loss of a gas is the relaxation frequency. The relaxation frequency is 
defined as that frequency at which the maximum vibrational absorption 
loss per unit wavelength occurs. In general, the relaxation frequency 
is a function of temperature, pressure, and humidity for a given gas. 

Air is assumed to be composed of nitrogen and oxygen, neglecting the 
absorption of the other components. Therefore, the following empirical 
expressions for the relaxation frequencies are used: 

^rl,n " (P/P^) (293.15/T)^/2^g + 35Qh exp|-6. 142 [(293 . 15/T) - l]|) (3) 

and 

^ = (p/Pj.)|24 + 44100h[(0.05 + h)/(0.391 + hj] I (4) 


In equations (3) and (4) , f„, ^ is the relaxation frequency of nitrogen 

in hertz, ^ is the relaxation frequency of oxygen in hertz, p is 

the ambient pressure, T is the ambient temperature in Kelvin, and h is 
the absolute humidity in percent mole fraction. Rewriting equations (3) 
and (4) in terms of dimensionless variables yields 

= p*/(T*)^/^^9.08 + 340. 65h exp |-6 . 178 [(T* ) - l] (5) 

and 

= P*|24 + 44100h[(0.05 + h)/(0.391 + hfj| (6) 

where T* = T/T^ and p* = p/p^* A graph of the relaxation frequency as 
a function of temperature and humidity at standard sea level pressure is 
presented in figure 1 for nitrogen. Relaxation frequency as a function of 
humidity at standard sea level pressure is presented in figure 2 for 
oxygen. It is interesting to note that the value of the oxygen relaxation 
frequency is always much greater than the nitrogen relaxation frequency 
and that the value of the relaxation frequency is highly dependent on the 
value of the absolute humidity. 
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The four terms on the right-hand side of equation (2) are now com- 
puted. The classical and rotational terms are combined and expressed as 
a function of temperature , pressure, and frequency as follows: 


“cl ^ “rot " (T/293.15)^/2 


Each of the two vibrational loss terms are written in the following form: 


a 


vib, i 



exp(-9j/T) (f/c) 2ffj,i^j^ 
jl - exp(-0j^/T)J ^ 


( 8 ) 


where 0^^ is the characteristic vibrational temperature and is the 

fractional molar concentration. In equation (8) , i is defined such that 
i = o for oxygen and i = n for nitrogen. Since the value of 0j^/T is 
large, the term in square brackets can be replaced by unity. Expressing 
equations (7) and (8) in terms of dimensionless variables, substituting 
the values for the physical constants given in table I, and using the 
relation 


c 




(9) 


yields the following: 


“cl “rot " (6.207 X 10-9) (f/c^) (T*)^/2(f/p*) 

“vib,o = (9.555 X 10-4) (f/Cj,) (T*)"^/2 
X exp[7.771(T* - DA*] 

and 


(10) 


( 11 ) 


“vib,n 


(1.683 X 10"4) (f/c^) (T*) 

X e.p[ll.633(T* - l)/T*]r2ff^^_„/(£2 . ^ J 


( 12 ) 


The total absorption coefficient is the sum of equations (10) , (11) , and 

(12) , The absorption coefficients in these equations are expressed in the 
units of nepers per meter. To convert to decibels per meter, it is neces- 
sary to multiply the right-hand side of each equation by 8.69. 
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Figure 3 shows a typical graph of the total absorption coefficient 
as a function of frequency with h=0.2, T*=1.0, and p* = 1.0. The 

three distinct regimes for the absorption coefficient are readily appar- 
ent from the figure. The first regime r where the frequency is less than 
the relaxation frequency of nitrogen, is dominated by the vibrational 
absorption of nitrogen. The second regime includes values of frequency 
between the nitrogen relaxation frequency and the oxygen relaxation fre- 
quency and is dominated by the oxygen vibrational absorption. The clas- 
sical and rotational losses dominate in the third regime for frequencies 
above the oxygen relaxation frequency. 

The data in figure 3 include a wide range of frequencies. For air- 
craft noise problems, the frequency range of interest is normally limited 
to less than 10 OQO Hz. To further demonstrate the properties of the 
absorption coefficient, the classical and rotational, nitrogen vibra- 
tional, oxygen vibrational, and total coefficients are plotted as func- 
tions of frequency in figures 4 to 7. The effect of changing relaxation 
frequency is shown by the lines of constant humidity in figures 5 to 7. 

All data are for standard sea level temperature and pressure, and all four 
figures are plotted to the same scale to allow direct comparisons between 
the figures. Finally, the relaxation frequencies which correspond to each 
constant humidity line are shown in figures 5 to 7. 

The total absorption-coefficient curves on figure 7 demonstrate how 
the characteristics of the absorption coefficient dramatically change as a 
result of changing relaxation frequency. As the values of the atmospheric 
properties change, the dominant absorption term that falls within the 
frequency range of interest changes. Comparison of figure 7 with figure 3 
can assist in identifying the dominant term. 

The atmospheric-absorption coefficient is dimensional. To express 
the coefficient in dimensionless form, a dimensionless absorption coeffi- 
cient p is defined as 


]i = acj./f 


(13) 


where a is the average absorption coefficient from the ground to a given 
height. The dimensionless absorption coefficient expresses the absorption 
loss in terms of nepers per unit wavelength under standard conditions. In 
general, y is a function of temperature, pressure, humidity, and fre- 
quency- For the Atmospheric Module used in ANOPP (ATM), the temperature, 
pressure, and humidity are a function of altitude only. The dimensionless 
altitude y, as discussed in the ATM, is defined as 


y = Mg^(H - 


(14) 


The Atmospheric-Properties Table from ATM gives temperature T*, 
pressure p , and humidity h as functions of y. Therefore, the 
absorption coefficient is expressed as a function of y and f. 
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The total change in sound intensity due to atmospheric absorption is 
an integral over the length of the conical ray tube from the source to 
the observer. An average absoirption coefficient is defined as this inte- 
gral divided by the length. Since the absorption coefficient is a func- 
tion of y and f only, the average dimensionless absorption coefficient 
from ground level to some value y is given by 


y(f,y) I (ctCj./f) dy (15) 

^ •'o 

Equation (15) is also valid for each component of the absorption coeffi- 
cient. Equations (10), (11), and (12) are then rewritten as 


V^cl + ^rot = 


/• y 

.207 X 10"^)- I (T*)^/Vp* dy 
^ •'O 


and 


/" V 2f f 

“vib.o - <5-555 X 2 ''I 2 

•^0 f + f„, 

X exp [ 7 . 771 (T* - 1)/T*] dy 


rl,o 


i-vib.n = I1-5S3 X lO-'l^ (T*)-5/2 

rl, n 


f 

« Jo 

X exp[ll.633(T* - 1)/T*] dy 


(16) 


(17) 


(18) 


The temperature T* , pressure p* , and relaxation frequencies are func- 
tions of y- The total average dimensionless absorption coefficient is 


y(f,y) + ^rot ^vib,o ^vib,n 

This module produces a table of y(f,y) for a range of values of the 
altitude H and frequency f for use by the Propagation Module. The 
Propagation Module corrects the sound intensity values for atmospheric- 
absorption effects. 
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REFERENCES 
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Proposed New Method for Computing Atmospheric Absorption Losses. 
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TABLE I.- STORED PRIMARY CONSTANTS 


Constant 

SI Units 

U.S. Customary 
Units 



340.294 m/s 

1116.45 ft/s 



9.806 65 m/s^ 

32.1741 ft/s^ 

M 

28.9644 

28.9644 

Pr 

1.225 kg/m^ 

0.002 377 slug/ft^ 

R 

8314.32 m^K-s^ 

49 718.96 ft^/°R-s^ 



288.15 K 

518.67°R 



0.781 

0.781 



0.209 

0.209 

1 

0n 

3352.0 K 

6033. 6°R 



2239.1 K 

4030. 38°R 

o 
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Relaxation Frequency, f^l ^ 





Relaxation Frequency 



Total Absorption Coefficient, a , nepers/meter 



Figure 3 .- Typical total adDSorption coefficient for standard 
sea level temperature and pressure. 
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Classical Plus Rotational Absorption Coefficient, nepers/meter 



Figure 4.- Classical plus rotational absorption coefficient for standard 
sea level temperature and pressure. 
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vibrational Absorption Coefficient, nepers/meter 



Figure 5.- Nitrogen vibrational absorption coefficient for standard 
sea level temperature and pressure. 
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Figure 6.- Oxygen vibrational absorption coefficient for standard 
sea level temperature and pressure. 
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Total Absorption Coefficient, «, nepers/meter 



Figure 7.- Total absorption coefficient for standard 
sea level temperature and pressure. 
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3.2 GROUND REFLECTION AND ATTENUATION MODULE 


INTRODUCTION 

The noise produced by aircraft is of most concern during take-off 
and landing, when the observer and the aircraft are both close to the 
ground. The ground reflection and attenuation effects are a significant 
influence on the propagation of sound waves under these conditions. The 
presence of the Earth's surface causes a change in the sound spectrum due 
to the reflection and attenuation of the sound waves and the creation of 
surface waves. The ground reflection and attenuation must be properly 
accounted for in order to accurately predict the aircraft noise at the 
observer. 

Sound waves which propagate nearly parallel to the Earth's surface 
may be attenuated by the absorption of the ground. The amount of atten- 
uation is greatly dependent on the surface characteristics. A surface is 
characterized by a complex acoustic impedance. The sound waves can also 
be reflected by the ground and can either enhance or diminish the sound 
intensity received directly, depending on the phase-angle shift. Finally, 
surface waves can be produced in addition to the normal free-field spheri- 
cal acoustic waves due to the ground influence. All of these effects must 
be accounted for to produce a complete model. 

The ground-effects model chosen for ANOPP (ref. 1) is the Chien- 
Soroka theory (ref. 2). This theory, with the Delany and Bazley impedance 
function (ref. 3) , assumes that the ground is a locally reacting uniform 
plane and that the source is a point source. This model has been vali- 
dated (ref. 4) with the data of Parkin and Scholes (refs. 5 and 6). 

The Ground Reflection and Attenuation Module computes a table of the 
ground-effects factor, which is a dimensionless coefficient incorporating 
the effects of reflection, attenuation, and surface waves. The ground- 
effects factor is tabulated in dimensionless form as a function of path- 
length difference, incidence angle, frequency, and source-to-observer 
distance. 


SYMBOLS 

a incoherence constant 

C coherence coefficient 

c speed of sound, m/s (ft/s) 

F spherical-wave shape function 

f frequency, Hz 
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G 

H 

h 

i 

K 

k 

% 

<p2> 

R 


ground-effects factor 
source altitude, m (ft) 
observer altitude, m (ft) 

unit imaginary number 

^l/(6Nu.) 
constant, 2 ^ 

wave number, 27Tf/c 

number of subbands per 1/3-octave band 
number of ground dips 

2 2 4 

mean-square pressure. Pa (lb /ft ) 


magnitude of complex spherical- wave reflection coefficient 
r distance, m (ft) 

U unit step function 

a argument of complex spherical-wave reflection coefficient 

r complex plane-wave reflection coefficient 

Ar path-length difference, m (ft) 

e constant, K - 1 

n = 27Tpf/a 

9 incidence angle, deg 

V complex specific ground admittance 

p density, kg/m^ (slugs/ft^) 

o 3 

a specific flow resistance of ground, kg/s-m^ (Ib/s-ft ) 

T = (kr2/2i) (cos 0 + V) 

Subscripts : 

c center band value 

ff free field 

gr ground effect 

Z lower limit 
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u 


upper limit 


1 direct 

2 reflected 


INPUT 


The basic independent variables for the model are path-length dif- 
ference, incidence angle, frequency, and source-to-observer distance. 

The range of path-length differences is computed within the program based 
on the user-specified number of ground dips. The ranges for other vari- 
ables are specified by upper and lower limits. The atmospheric properties 
of density and speed of sound and the specific flow resistance of the 
ground are also required. These values are assumed to be constant 
throughout the model. The number of 1/3-octave siabband intervals adjusts 
the predicted effect for bandwidth. The incoherence constant is an 
empirical quantity which limits cancellation effects. The range and 
default values of each input are given in table I. 


a 

c 


% 

Nd 

P 


a 


incoherence coefficient 

speed of sound at the observer, m/s (ft/s) 
frequency lower and upper limits, Hz 
number of subbands per 1/3-octave band 
number of ground dips to be included 

source-to-observer distance lower and upper limits, m (ft) 

incidence angle lower and upper limits, deg 

air density at the observer, kg/m^ (slugs/ft^) 

specific flow resistance of the ground, kg/s-m^ (slugs/s-f t^) 


OUTPUT 

The module produces a tcd^le of the ground-effects factor as a func- 
tion of four dimensionless parameters: path-length difference, cosine of 

incidence angle, frequency, and source- image- to-observer distance. 


Ground Effects Table 
kAr path-length difference 

cos 9 cosine of incidence angle 
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n 


= 27 Tpf/a 


kr 2 image distance 

G(kAr,cos 9,r)rkr2) ground-effects factor 

METHOD 

The ground-effects geometry is shown in figure 1. A source is 
located at an altitude H over a ground plane. Sound arrives at the 
receiver at a height h from the direct path and from a reflected 

path r 2 ^ which appears to the observer to be from an image source. The 
incidence angle of the reflected wave is 9. The path-length difference 
Ar = T 2 “ is the most significant parameter of ground effects. As 

shown in figure 2 , the path- length difference can be approximated in 
terms of the observer height and the incidence angle as 


Ar - 2h cos 9 


( 1 ) 


The Chien-Soroka theory is derived from a solution to the wave equa- 
tion in the half space of figure 1. The derivation of the theory is pre- 
sented in references 1 and 2. The resulting expression for the mean- 
square pressure with ground effect is 

<p^>gr = <P^>ff t + 2RC COS (a + kAr)] (2) 

where is the free-field mean-square pressure, C is the coher- 

ence coefficient, k is the wave number, R is the magnitude, and a 
is the argument of the complex spherical-wave reflection coefficient. The 
term in brackets in equation (2) is referred to as the ground-effects 
factor G and is defined as follows: 


G - <P^>grAp^>ff 


(3) 


The coherence coefficient C is the fraction of the initial acoustic 
energy in which phase relation is maintained throughout the propagation 
process. A reasonable approximation for the coherence coefficient is made 
by assuming a Gaussian distribution of the form 


C = exp[^{akAr)^ 


( 4 ) 


where a is the incoherence constant and exp denotes the exponential 
function. The incoherence constant is normally given a value of 0.01, 
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which corresponds to a value of C of 0.37 at a Ar value of 16 wave- 
lengths. After substitution of equation (4), the ground-effects factor 
becomes the following: 


G 


1 + + 2R exp (akAr) ^ cos (a + kAr) 


( 5 ) 


The Chien-Soroka theory (ref. 2) shows that the complex spherical- 
wave reflection coefficient can be expressed as follows: 


1 nt 

Re-" = r + (1 - DF(t) 


( 6 ) 


where F is the complex plane-wave reflection coefficient given as 


r = (cos 6 - v)/(cos 0 + V) 


(7) 


and F(t) accounts for the spherical-wave shape. In this equation V is 
the complex specific ground admittance. The function F(T) in equa- 
tion (6) is 


F(t) = 1 - JVt W(iT) 


(8) 


where 


T = (kr2/2i) (cos 6 + V) 


(9) 


and W is the following complex error function: 


W(z) 



(lm(z) > 0) 


( 10 ) 


For any value of [t| > 10, an asymptotic approximation for the complex 
error function is used. This allows F(t) to be expressed as the 
following : 


F(t) 


-2 \/Tru[[-Re(T)] Te^^ + 


2 2 
(2t2)^ 


( 11 ) 
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where U is the unit step function defined as follows: 


U(S) =1 
U(S) = 1/2 
U(S) =0 


(S > 0) 

(S = 0) 

(S < 0) J 


(12) 


The remaining parameter to be determined is the complex specific 
ground admittance V. Delaney and Bazley (ref. 3) developed the following 
empirical equation for V: 

V = [l + {6.86n)”°'^^ + i(4.36n)"°'^^[l ^ (13) 

where the dimensionless frequency n is 


ri = 2irpf/a 


(14) 


A graph of the ground admittance is shown in figure 3. 

For an acoustically hard surface (n = 0) , the theory is greatly sim- 
plified since n=0, T=l, R=l, and a = 0. The expression for the 
ground-effects factor reduces to 


G = 2 + 2 exp|j- (akAr) cos (kAr) 


(15) 


which is a function of kAr only. 

Actual noise predictions are given for finite- frequency bandwidths. 
In order to approximate the effect of finite bandwidth, the ground- 
effects factor is integrated over the band, but only the variation of the 
term cos (a + kAr) is included as a variable in the integration. It is 
assumed that the variation of the other terms in the width of 1 band is 
small. This approximation is acceptable for 1/3-octave or narrower 
bands (ref. 4) . 

For purposes of analyzing attenuation and reflection effects, the 
standard 1/3-octave bands of sound are subdivided into Nj^ subbands, 
where Nj^ is an odd number. Using an odd number gives a center subband 
which has a center frequency equal to the center frequency of the 
original 1/ 3-octave band. 

. . ^ . . ^l/(3Nh) ... .. ..V, 

The ratio of the subband limit frequencies is 2 , so that the 

integral for averaging the cosine term over the sxabband is the 
following : 
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cos (a + kAr) df 


(16) 


<cos 


1 

(a + kAr)> = I 

Jf/K 


1/(6N^) 

where K = 2 ^ , Af = (K - and is the center frequency 

of the subband. Assuming that a remains constant throughout the sub- 
band, an approximation for equation (16) is expressed as 

sin (ek Ar) 

<cos (a + kAr)^ = cos (a + k^Ar) - v (17) 


where k^ is the subband center wave number and e = K - 1. This pro- 
duces the following final expression for the ground-effects factor: 


(18) 


where it is understood that k refers to the subband center wave number. 
For an acoustically hard surface, the averaged expression for the final 
band is 


G = 2 + 2|exp^- (akAr) cos (kAr) 


(19) 


A plot of G as a function of kAr is presented in figure 4 for the case 
of an acoustically hard surface. 

It is apparent from figure 4 that the ground-effects factor has a 
series of nodes occurring at values of kAr of (2n - 1)tt, where 
n=l, 2, ..., n+1. These nodes are caused by the cancellation of 

sound intensity between the direct and reflected waves and are commonly 
referred to as ground dips. It is imperative that the value of G at 
the minimum of each ground dip as well as at intermediate points be 
evaluated to adequately define the function; however, after the fifth dip 
(Njj = 5) the variation in G for increasing values of kAr is so small 
that it can be assumed to be constant. If this condition is not met, the 
user can provide a different value for N^. In addition, the frequency, 
source- to-observer distance, and incidence-angle range of interest are 
user inputs. These three variable ranges are divided into equally spaced 
intervals. 

For most conditions of interest for which the noise-generation models 
are valid, the source and observer separation is large. This means that 
h/r is much less than unity. Under this condition it can be assumed 
that T 2 - and Ar - 2h cos 0 2 hH/r 2 ^. 
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The output of the module is a table of the ground-ef fects-factor 
values as functions of four dimensionless variables: kAr, cos 9/ T\, 

and kr 2 « For a hard surface, the last three variables do not affect the 
ground-effects factor. In this case, the module output is a four- 
dimensional table with only one entry in each of the dimensions cos 0, 
r)/ and kr 2 * The interpolation' logic produces the same result for this 
table as it would for a one-dimensional table in kAr. 

In reference 4, the ground effects predicted by this method have 
been compared to the ground-to-ground propagation data from references 5 
and 6. The predictions and data are in fair to good agreement in this 
case. Recent comparisons of this method with T-38A flyover data given 
in reference 7 and Boeing 747 flyover data given in reference 8 show 
that the predicted frequency of maximum ground attenuation agrees very 
well with the measured frequency. The amplitudes of the predicted 
attenuation are generally smaller than the measured amplitudes, indi- 
cating that the Chien-Soroka method may be a conservative estimate of 
the effects of ground reflection and attenuation. 
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TABLE I.- RECOMMENDED RANGES FOR INPUT PARAMETERS 


1 

Input 

— 

Minimum 

Default 

Maximum 

^ • • • • • 

• 

. 

. 

0.001 

0.010 

0.100 

c, m/s . . . 



• 

300 

340.294 

400 

f£' • • • 




13 

50 


Hz . . . 





2000 

4000 

Nfa 




1 

5 

9 

Nd 




2 

5 

10 

r£' m ... 




10 

10 


r^, km . . . 





10 

10 

0£, deg . . 




0 

0 


6^, deg . . 





89 

89 

p, kg/m^ . . 




1.0 

1.225 

1.5 

3 

a, kg/s-m 

• 

• 


1.0 X 10^ 

2.5 X 10^ 

5.0 X 10^ 
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Figure 1.- Ground-effects geometry. 
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Ground Effects Factor 



Figure 4.- Ground-effects factor for acoustically hard surface. 
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4. SOURCE NOISE PARAMETERS 




4,1 FAN NOISE PARAMETERS MODULE 


INTRODUCTION 

Fan noise is a significant part of the total noise produced by turbo- 
jet and turbofan engines. Fan noise prediction methods are included in 
ANOPP, The purpose of the Fan Noise Parameters Module is to generate the 
physical parameters required by Heidmann's method (ref. 1) for fan noise 
prediction for turbojet and turbofan engines. 


The fan entrance and exit flow states are provided by the user. 
These engine-state data are first converted to a function of time, using 
data from the flight-path characteristics. The required fan parameters 
for fan noise prediction are then computed. 


A 

c 

D 

M, 

m 

N 

n 

P 

R 

T 

t 

Y 

7T 


SYMBOLS 

area, m^ (ft^) 

speed of sound, m/s (ft/s) 

fan diameter, m (ft) 

aircraft Mach number 

mass flow rate, kg/s (slugs/s) 

rotational speed, Hz 

number of source times 

pressure. Pa (Ib/ft^) 

dry-air gas constant, m^K-s^ (ft^/^R-s^) 
temperature, K (^R) 
time, s 

ratio of specific heats 
power setting 


P 


density, kg/m^ (slugs/ft^) 



Subscripts : 


e 

engine 

i 

entrance 

j 

exit 

t 

total 

oo 

ambient 


Superscript; 

* dimensionless quantity 


INPUT 

This module provides the parameters for a typical axial-flow fan as 
shown in figure 1. The entrance and exit flow states for the fan are 
required from the user. The engine power setting, aircraft Mach number, 
and ambient density and speed of sound are provided by the engine variable 
table. 


7T 

M 




N (TT,M^) 


Fan Entrance Flow State 
engine power setting 
aircraft Mach number 
mass flow rate, re 
rotational speed, re ^RT^^D 


T^ ^(tt,M^) total temperature, re T^ 


7T 


M 


m. (TT,MJ 


T^^ .(7T,MJ 


Fan 

engine power setting 
aircraft Mach number 
mass flow rate, re 
total temperature , 


Exit Flow State 


p A p A/RT 
re T 

oo 


Engine Variable Table 


t source time, s 

M^(t) aircraft Mach number 


4.1-2 



7T(t) 


engine power setting 


c^(t) ambient speed of sound, m/s (ft/s) 

Poo(t) ambient density, kg/m^ (slugs/ft^) 


OUTPUT 


The outputs 
execution of the 


to this module are the physical parameters required for 
fan noise modules as a function of source time. 


n 

t 

in. (t) 

1 

m*(t) 

N*(t) 

At* ( t) 


Fan Noise Parameters 
number of source time values 
source time, s 

entrance mass flow rate, re p c A 

' ^00 oo e 

exit mass flow rate, re p c A 

QO 00 e 

rotational speed, re c^D 

total temperature rise across fan, re T 

oo 


Ambient Conditions 

c^(t) ambient speed of sound, m/s (ft/s) 

M^(t) aircraft Mach number 

p^(t) cunbient density, kg/m^ (slugs/ft^) 


METHOD 

The fan entrance and exit flow states are expressed as a function of 
the engine power setting 7T and the aircraft Mach number M , These 
tables must be provided directly by the user. These data are converted 
to a function of source time. The engine variable table gives M and 7T 
as a function of source time as provided by the Flight Dynamics Mo^dule or 
the user. 

The fan entrance and exit mass flow rate and the fan rotational speed 
must be converted to the referred variables used by the Fan Noise Module 
by the relations 


m (t) 


m 




( 1 ) 
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and 


N*(t) 


4r 


N 




( 2 ) 


where the ambient value of the ratio of specific heats Y is 1.4. 
Finally, the temperature rise across the fan At* is the difference in 
total temperature. 


. * 

AT 




(3) 


REFERENCE 

1. Heidmann, M. F. : Interim Prediction Method for Fan and Compressor 

Source Noise. NASA TM X-71763, 1975. 
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4.2 CORE NOISE PARAMETERS MODULE 


INTRODUCTION 

Core noise is a significant part of the total noise produced by 
turboprop, turbo fan, and turbojet engines. Core noise prediction methods 
are included in ANOPP, The purpose of the Core Noise Parameters Module 
is to generate the physical parameters required for a core noise predic- 
tion module. 

The core entrance and exit flow states are provided by the user. 
These engine-state data are first converted to a function of time, using 
data from the flight-path characteristics. The required core parameters 
for core noise prediction are then computed. 


SYMBOLS 

A area, m^ (ft^) 

c ambient speed of sound, m/s (ft/s) 

M aircraft Mach number 

oo 

m mass flow rate, kg/s (slugs/s) 

n number of source times 

p pressure. Pa (Ib/ft^) 

R dry-air gas constant, m^/K-s^ (ft^/^R-s^) 

T temperature, K (^R) 

t time, s 

7T engine power setting 

3 3 

p density, kg/m (slugs/ft ) 

Subscripts : 
e engine 

i entrance 

j exit 


4.2-1 



t 


total 


00 ambient 

Superscript: 

* dimensionless quantity 


INPUT 

The entrance and exit combustor flow states are required from the 
user. The engine power settings and aircraft Mach numbers in the engine 
variable table are provided by the Flight Dynamics Module or the user. 


Core Entrance Flow State 


7T 

engine power setting 


M 

oo 

aircraft Mach number 


m* ) 

1 oo 

mass flow rate, re 

A p /JrT 

e'^oo /V oo 


total pressure, re 

Poo 


total temperature. 

re 


Core Exit Flow State 
TT engine power setting 

M aircraft Mach number 

00 

T* total temperature, re 

Engine Variable Table 

source time, s 
aircraft Mach number 
engine power setting 
cunbient speed of sound, m/s (ft/s) 
ambient density, kg/m (slugs/ft ) 

OUTPUT 

The outputs to this module are the physical parameters required for 
the execution of the core noise modules as a function of source time. 


M^(t) 

00 


TT(t) 

Cco(t) 

P^(t) 
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Core Noise Parameters 


n number of source time values 

t source time r s 

m*(t) combustor entrance mass flow rate, re p c A 

p* (t) combustor entrance total pressure, re p^ 

T* (t) combustor entrance total temperature, re T 

1 

T*(t) combustor exit total temperature, re T 

D 

Ambient Conditions 

c (t) ambient speed of sound, m/s (ft/s) 

oo 

M (t) aircraft Mach number 

oo 

3 3 

p (t) cunbient density, kg/m (slugs/ft ) 

oo 


METHOD 

A schematic diagram of a typical combustor depicting the entrance 
and exit flow states is shown in figure 1- These flow states are a func- 
tion of the engine power setting tt and the aircraft Mach number M^. 
These data are converted to a function of source time by interpolating 
the core flow state tables with respect to M^(t) and tt (t) values from 
the engine variable table. In addition, the combustor entrance mass flow 
rate must be converted to the referred variables used by the Combustion 
Noise Module by the relation 


(t) 



( 1 ) 


where the ambient value of the ratio of specific heats 


Y 


is 1.4. 
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4.3 TURBINE NOISE PARAMETERS MODULE 


INTRODUCTION 

Turbine noise is a significant part of the total noise produced by 
turbofan and turbojet engines during low-power operations. Turbine noise 
modules are included in ANOPP. The purpose of the Turbine Noise Param- 
eters Module is to generate the physical parameters required for turbine 
noise prediction for turbojet and turbofan engines. 

The turbine entrance and exit flow states are provided by the user. 
These engine-state data are first converted to a function of time, using 
data from the flight-path characteristics. The required turbine param- 
eters for turbine noise prediction are then computed. 

SYMBOLS 

A area, m^ (ft^) 

c speed of sound, m/s (ft/s) 

D turbine rotor diameter, m (ft) 

f fuel-to-air ratio 

h* specific enthalpy, re Rt 

00 

h absolute humidity, percent mole fraction 

a 

M Mach number 

M aircraft Mach number 

oo 

m mass flow rate, kg/s (slugs/s) 

N rotational speed, Hz 

n number of source times 

p pressure. Pa (Ib/ft ) 

R dry-air gas constant, m^/K-s^ (ft^/^R-s^) 

R gas constant, m^K-s^ (ft^/^R-s^) 

T temperature, K (°R) 

t time , s 


4.3-1 



Y ratio of specific heats 

TT engine power setting 

p densityr kg/m'^ (slugs/ft ) 

({)* specific entropy function, re R 

Subscripts : 
e engine 

i entrance 

j exit 

s static 

t total 

°° ambient 

Superscript: 

* dimensionless quantity 


INPUT 

This module provides the parameters for a typical axial-flow turbine 
as shown in figure 1. The entrance and exit flow states for the turbine 
are required from the user for predicting the turbine noise. The engine 
power settings, aircraft Mach numbers, and ambient densities are provided 
by the engine variable tables. 


M 


m. (7T,M ) 

1 ' * co' 


N (TT,M ) 
00 

* 


Turbine Entrance Flow State 
engine power setting 
aircraft Mach number 

mass flow rate, re o/\/^ 
rotational speed, re 


Tt ^(TT,M^) total temperature, re T^ 


Turbine Exit Flow State 
TT engine power setting 

M_ aircraft Mach number 

oo 
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J °° 

f (TT,M^) 
00 


m. (TT»M ) 

J °° 


* 



j 



(7T,M ) 
' co' 

(’^,M ) 

00 


turbine exit area^ re 
fuel-to-air ratio 

exit mass flow rate, re /V^^oo 

exit total pressure, re p 

'*■ oo 

exit total temperature, re T 

00 


Engine Variable Table 

t source time, s 

M^(t) aircraft Mach number 

TT(t) engine power setting 

c^(t) ambient speed of sound, m/s (ft/s) 

p^(t) ambient density, kg/m^ (slugs/ft^) 

hg(t) absolute humidity, percent mole fraction 


OUTPUT 

The outputs of this module are the physical parameters required for 
the execution of the turbine noise modules as a function of source time. 


Turbine Noise Parameters 


n 

number of source time values 

t 

source time, s 


f(t) 

fuel-to-air ratio 


m* (t) 

mass flow rate, re p c A 

oo 00 e 


* . . 

N (t) 

rotational speed, re c 



entrance total temperature, 

re T 
00 


exit static temperature, re 

T 

oo 


Ambient Conditions 


ambient speed of sound, m/s 

(ft/s) 

ha(t) 

absolute humidity, percent mole fraction 
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M (t) aircraft Mach number 

oo * 

p (t) ambient density, kg/m (slugs/ft^) 

00 

METHOD 

The turbine entrance and exit flow states are expressed as a function 
of the engine power setting tt and the aircraft Mach number M^. These 
tables are provided directly by the user. These data are converted to a 
function of source time. The engine variable table gives M^ and 7T as 
a function of source time as provided by the Flight Dynamics Module or 
the user. 

The turbine rotational speed and mass flow rate must be converted to 
the referred variables used by the turbine noise modules by the relation 

N*(t) = ^ (1) 

and 

m*(t) = 5!— (2) 

\n 

where the ambient value of the ratio of specific heats y is 1.4. The 

turbine exit static temperature T . is required by the turbine noise 

s , 3 

modules for the computation of the ideal work extraction. The exit static 
temperature can be computed assuming either constant specific heats or 
variable specific heats. 


Constant Specific Heats 

The ratio of specific heats at the turbine exit y^ and the dimen- 
sionless gas constant R* = R/R are computed from the total temperature 

T* , fuel-to-air ratio f, and absolute humidity h , using the 
t, j 

appropriate thermodynamic utility. The turbine exit Mach number M^ 
can be computed by the relation 




1 ^ 
2 Y ^-1 


( 3 ) 
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as discussed in Thermodynamic Utilities. Then, the turbine exit static 
temperature T . is 
s,i 



Variable Specific Heats 

As discussed in Thermodynamic Utilities, the static temperature T* 
can be computed from simultaneous solution of 


and 


_ 


tr j 


= e 




( 5 ) 


,/rt* . 

3 V t,T 

TV* * 

A , p^ ■ 

: to 


111 


so ^to 


- <) 


( 6 ) 


where <() xs the dimensionless entropy function and h* is the dimen- 
sionless enthalpy. 
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4.4 JET NOISE PARAMETERS MODULE 


INTRODUCTION 

Exhaust jet noise is a significant part of the total noise produced 
by turbojet and turbofan engines. Several jet noise prediction methods 
are included in ANOPP. The purpose of the Jet Noise Parameters Module is 
to generate the physical parameters required for the exhaust jet noise 
prediction modules for turbojet and turbofan engines. 

The engine nozzle exit flow state is provided by the user. These 
engine-state data are first converted to a function of time, using data 
from the flight path characteristics. The required jet parameters for jet 
noise prediction are then computed. 


SYMBOLS 

A area, m^ (ft^) 

c speed of sound, m/s (ft/s) 

equivalent circular nozzle diameter, m (ft) 
hydraulic diameter, m (ft) 

Dp plug diameter, m (ft) 

f fuel-to-air ratio 

absolute humidity, percent mole fraction 
M Mach number 

M^ aircraft Mach number 

m mass flow rate, kg/s (slugs/s) 

n number of source times 

p pressure. Pa (Ib/ft^) 

R dry-air gas constant, m^K-s^ (ftV°R-s^) 

R gas constant, m^K-s^ (ft^/^R-s^) 

T temperature, K (®R) 

t time , s 
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V 


velocity, m/s (ft/s) 


Y ratio of specific heats 

7T engine power setting 

p density, kg/m (slugs/ft ) 

(j)* specific entropy function, re R 

Subscripts : 
a aircraft 

e engine 

fe fully expanded 

P plug 

s static 

t total 

1 primary strecim 

2 secondary streaon 

a> ambient 

Superscript; 

* dimensionless quantity 


INPUT 

This module provides the parameters which are used for the predic- 
tion of jet mixing and shock cell noise. The primary nozzle flow state 
is required for single-stream jets such as on turbojet engines. Both the 
p]fxmary and secondary nozzle flow states are required for dual— stream jets 
such as on turbofan engines. These jets and their associated state 
variables are illustrated in figure 1. The areas are included as state 
variables since the nozzle sizes may be functions of engine power setting. 
In addition to the areas and A 2 , the primary nozzle may have a plug 

with area A^, which must be specified, and the engine reference area A^ 
must be given since input data will usually be specified relative to it. 

It is assumed that the outer diameter of the inner nozzle is equal to the 
inner diameter of the outer nozzle as shown in figure 1. Local atmo- 
spheric conditions - density, speed of sound, and humidity - are used in 
evaluation of the static flow variables; however, humidity may be 
neglected without significant error if it is not available. The local 
density and speed of sound are used in computing decibel levels; however, 
they are not needed if only normalized mean-squared pressure is to be 
computed or if the SPL is corrected to standard conditions. Table I 
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gives the recommended ranges and the default values for the input 
parameters . 

2 2 

engine reference area^ m (ft ) 

A* primary nozzle plug area, re A 

p 

Primary-Nozzle Flow State 

7T engine power setting 

M aircraft Mach number 

00 

A*(7T) nozzle area, re A^ 

fuel-to-air ratio 

m*(7T,M^) mass flow rate, re A^p^^^RT^ 

it 

Pt total pressure, re p^ 

T* total temperature, re 

Secondary-Nozzle Flow State (Optional) 

7T engine power setting 

M aircraft Mach number 

00 

it 

A (7T) nozzle area, re A 

f 2 ('n',M^) fuel-to-air ratio 

m*(7T,M^) mass flow rate, re A^p^^,JrT^ 

p* total pressure, re p^ 

* 

T. ^('IT,M ) total temperature, re T 
t , 2 °° 

Engine Variable Table 

t source time, s 

M (t) aircraft Mach number 

00 

TT (t) engine power setting 

c^(t) cunbient speed of sound, m/s (ft/s) 

o o 

p^(t) ambient density, kg/m (slugs/ft ) 
h^(t) absolute humidity, percent mole fraction 
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OUTPUT 


The outputs to this module are the physical parcimeters required for 
execution of the exhaust jet noise modules as a function of source time. 
Unless otherwise stated, all paraimeters are computed for a hypothetical 
fully expanded jet which has a static pressure equal to ambient pressure. 


T*(t) 

V*(t) 


Pl(t) 

Yi(t) 


Primary Jet Parameters 
number of source time values 
source time , s 
primary jet area, re 

actual primary jet equivalent diameter, re yj 

actual primary jet hydraulic diameter, re yj A^ 

primary jet Mach number 

primary jet total temperature, re T^ 

primary jet velocity, re c^ 

primary jet density, re 

ratio of specific heats for primary jet 


°h, 2 ^^> 

M2<t) 

T;(t) 

v;(t) 

p*(t) 


Y2(t) 


Secondary Jet Parameters (Optional) 
secondary jet area, re A^ 

actual secondary jet equivalent diameter, re 

actual secondary jet hydraulic diameter, re \J ^ 

secondary jet Mach number 

secondary jet total temperature, re T^ 

secondary jet velocity, re c^ 

secondary jet density, re 

ratio of specific heats for secondary jet 
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Ambient Conditions 


c^(t) ambient speed of sound, m/s (ft/s) 

M- (t) aircraft Mach number 

a 

Q 'i 

p^(t) ambient density, kg/m (slugs/ft ) 


METHOD 

The user has two options for the method of this module. The first 
is the previously ■ stated option of computing data for either a single- 
stream or dual-stream jet. The second is the selection of constant ratio 
of specific heats or variable ratio of specific heats for calculation of 
static state variables. 

The input nozzle flow states as shown in figure 1 are expressed as a 
function of the engine power setting 7T and the aircraft Mach number M^. 
These tables must be provided directly by the user. These data are con- 
verted to a function of source time. The engine variable table gives M^ 
and 7T as a function of source time as provided by the Flight Dynamics 
Module or the user. 


Single-Stream Jet 

The ratio of specific heats for the primary jet is computed 

from the total temperature T , the fuel-to-air ratio f , and the 

t , 1 

absolute humidity h^, as discussed in the Thermodynamic Utilities portion 
of this manual. 

Constant specific heats .- The primary jet velocity and jet density 
may be estimated assuming constant specific heats. The jet Mach num- 
ber is computed from the jet total pressure and ratio of specific 

heats by the relation 




P 


s,l 






( 1 ) 


where p* = 1 for a fully expanded jet. 
s / 1 

The primary jet static temperature is computed from the jet Mach 
number and specific heat ratio as 



(2) 
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as discussed in Thermodynamics Utilities. The jet density is then 


* , . 





(3) 


where R* = R/R is the primary jet gas constant. Finally, the jet 
velocity is given by 


Y R*T* . 
s s, 1 


V* = mJ 


.1/2 


Y„ 


(4) 


where the ambient ratio of specific heats is 1.4. 

Variable specific heats .- As discussed in Thermodynamic Utilities, 
the primary jet static temperature is evaluated from the relation 


^s,l ^ 

* 

Pt,l 


where <})* is the dimensionless entropy function and p* = 1 for a 

5 / J- 

fully expanded jet. The jet Mach number is then given by 


M 


1 


I -rt* * • * 

^s,l "'1 


* * 

^l^s,! 


(6) 


where Y is the ratio of specific heats evaluated at T* . Finally, 
the jet density Pj^(t) and the jet velocity are computed from 

equations (3) and (4), respectively. 


Jet geometry .- The area of the fully expanded primary jet is 


A 


* 

fe,l 


** / *TT* 

mi/PiVi 


(7) 


The equivalent and hydraulic diameters, and are based on actual 

geometric values for the nozzles instead of the fully expanded area, which 
may occur slightly away from the nozzle exit plane. The equivalent 
diameter is 


“^1 ' \K/^ 
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and the hydraulic diameter is 



Secondary Jet 

Computations for the secondary stream are identical to those for the 
primary, except that the secondary plug diameter is taken to be the first 
outer diameter of the inner nozzle. Therefore, the hydraulic diameter of 
the secondary jet is 
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TABLE I.- RANGE AND DEFAULT VALUES OF INPUT PARAMETERS 


Input 
parame ter 

Minimum 

Default 

Maximum 

^ 2 
Ag /HI 

0.01 

it/ 4 

10 

A* 

0 

0 

1 

P ’ 








(a) Single-streeun circular nozzle. 



SECONDARY- NOZZLE 
FLOW STATE 
A^, m^, Pg. T^ 


Figure 1.- Schematic diagram of a circular and coannular 
nozzle depicting the appropriate flow states. 
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4.5 AIRFRAME NOISE PARAMETERS MODULE 


INTRODUCTION 

Airfraume noise is a significant part of the total noise of the air- 
craft during low-power operations. Airframe noise modules are included 
in ANOPP. The Airframe Noise Parameters Module generates the three time 
parameters needed for the execution of the airframe modules. The physical 
values output by the module are identical to the input; however, the input 
is in the form of a data member in ANOPP, whereas the prediction modules 
must have their input as parameters . 

SYMBOLS 

ambient speed of sound, m/s (ft/s) 
landing-gear position 
aircraft Mach number 
number of source times 
time, s 

flap setting, deg 

aimbient dynamic viscosity, kg/m-s (slugs/ft-s) 
ambient density, kg/m (slugs/ft ) 

INPUT 

Input to this module is the engine variable table as provided by the 
Flight Dynamics Module or the user. 

Engine Variable Table 

t source time , s 

I^^(t) landing-gear position 

M (t) aircraft Mach number 

oo 

6^(t) flap setting, deg 
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ambient speed of sound, m/s (ft/s) 


Poo(t) ambient density, kg/m^ (slugs/ft^) 

y^(t) ambient dynamic viscosity, kg/m-s (slugs/ft-s) 

OUTPUT 

The outputs to this module are the three time parameters required to 
execute the airframe noise modules as a function of source time. 

Airframe Noise Parameters 
number of source time values 
source time, s 
landing-gear position 
aircraft Mach number 
flap setting, deg 

Ambient Conditions 
ambient speed of sound, m/s (ft/s) 
ambient density, kg/m (slugs/ft ) 
ambient dynamic viscosity, kg/m-s (slugs/ft-s) 

METHOD 

The engine variable table is read and the appropriate output 
parameters are generated. 


Pco(t) 

Poo(t) 
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5. PROPAGATION 




5.1 PROPAGATION MODULE 


INTRODUCTION 

The Propagation Module (PRO) takes noise data which are generated by 
the noise modules in the source frame of reference and applies all of the 
appropriate computations to transfer them to the observer frame of refer- 
ence. If a sum of two or more noise sources is desired in the same 
coordinate system, the sum is performed first. The required computations 
for propagation include the effects of spherical spreading, atmospheric 
absorption, ground reflection and attenuation, and retarded time. The 

for some of the various propagation effects have been previously pre- 
pared by the Geometry Module (GEO) , the Atmospheric Absorption Module 
(ABS) , and the Atmospheric Module (ATM) . 

There are six steps necessary to propagate the sound to the observer. 
They are performed in the following order: 

1. Interpolate the input source noise data as a function of emission 

time and observer directivity angles 

2. Apply spherical spreading and characteristic impedance change 

effects 

3. Divide the 1/3-octave frequency bands into subbands for more 

accurate absorption and ground effects modeling (ref. 1) 

4. Apply the atmospheric absorption effect 

5. Apply the ground reflection and attenuation effect 

6. Combine frequency subbands into 1/3-octave bands 

The resulting output table is the mean— square acoustic pressure at the 
observer position as a function of frequency, reception time, and 
observer. 


SYMBOLS 

A subband amplitude adjusting factor 

c speed of sound, m/s (ft/s) 

f frequency, Hz 

Uj. acceleration due to gravity, m/s^ (ft/s^) 

H altitude, m (ft) 
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h 

k 

M 

o 

<P^> 

R 

r 

Ar 

SPL 

T 

t 

Ur V 
W 

y 

Y 

n 

0 

9 

M 

P 

a 

4) 

n 


observer height, m (ft) 

wave number, 27Tf/c^ 

molecular weight of air 

number of sxibbands per 1/3-octave band 

observer index 

2 2 4 

mean-square acoustic pressure. Pa (lb /ft ) 

universal gas constant, m^/K-s^ (ft^/^R-s^) 

distance, m (ft) 

path-length difference, m (ft) 

sound pressure level, dB 

temperature, K (^R) 

time, s 

spectrum slopes 

ratio of subband center frequencies 

dimensionless altitude, Mg^(H - h^)^RT^ 

elevation angle, deg 

= 27Tp^f/a 
a 

incidence angle, deg 
polar directivity angle, deg 

average absorption coefficient, nepers per wavelength 
density, kg/m^ (slugs/ft^) 

3 3 

specific flow resistance of the ground, kg/s-m (slugs/s-ft ) 
azimuthal directivity angle, deg 
solid angle for a ray cone, sr 


Subscripts : 
a ambient 


e 


emission 


o 


observer 
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r 


standard sea-level reference value 


s source 

Superscript: 

♦ dimensionless quantity 


INPUT 

The input to this module consists of one or more noise data tables 
computed for the saune source coordinate system on the aircraft. When two 
or more source noise data tables are input, they are sximmed prior to being 
propagated to the observer. Additional tables are required which incor- 
porate the various propagation effects. The range and default values of 
the required input parameters are given in tcible I. 

number of subbands per 1/3-octave band 

source radius, m (ft) 

specific flow resistance of the ground, kg/s-m^ (slugs/s-f t^) 
Source Noise Data Table 

f frequency, Hz 

0 polar directivity angle, deg 

(}) azimuthal directivity angle, deg 

tg emission time, s 

9 * 2 4 

<p^ (f ,0,<|),t^)^ mean- square acoustic pressure, re 


Geometry Table 


t 

reception time, s 

o 

observer index 

r (t,o) 

distance, m (ft) 

ft 

(P 

ft 

0 

emission time, s 

6(t,o) 

polar directivity angle, deg 

<|> (t,o) 

azimuthal directivity angle, deg 

Y(t,o) 

elevation angle, deg 

h(o) 

observer height, m (ft) 


^s 

a 
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Atmospheric Properties Table 


y altitude, Mgj.(H - Hj_)/rTj. 

c* (y) speed of sound, re c^. 

pc* (y) characteristic impedance, re 

p* (y) density, re 

Absorption Coefficient Table 
frequency, Hz 
altitude, Mg^(H - h^)^RT^ 

average absorption coefficient, nepers per wavelength 

OUTPUT 

This module produces a table of mean^square acoustic pressures at 
the observer as a function of frequency, time, and observer. 

Received Noise Data Table 

f frequency, Hz 

t reception time, s 

o observer index 

c*(o) speed of sound at the observer, re c^ 

a 

p* (o) air density at the observer, re 

\ 2 4 

<p^(f,t,o)> mean-square acoustic pressure, re p^c^ 

METHOD 

Noise Data Interpolation 

If two or more noise data tables are input, the tables are summed to 
produce one data table. Since the input noise data are expressed in 
terms of mean-square pressure, the noise data summation is a simple 
element-by-element addition. This noise data tcible is a function of 
source time values and constant polar and azimuthal directivity angle 
increments (pseudo observers) . Before these data can be propagated to 
the observer, they must be interpolated for the emission time values 
corresponding to reception time, polar directivity angles, and azimuthal 
directivity angles to the actual observers. This is accomplished by 


f 

Y 

U(f,Y) 
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using the polar directivity angle, azimuthal directivity angle, and 
emission time values from the geometry table. That is, using 0(t,o), 

()>(t,o), and tg(t,o), the input table <p^ (f ,0,(|),tg)> is converted 

to <p^(f,t,o)> . Figure 1 is a graphic representation of the time 
scale mapping. 


Spherical Spreading and Characteristic Impedance Effects 
2 ^ 

The table ^p (f,t,o)^ is now processed on a column-by-column 
basis throughout the remainder of the module for each value of received 
time and observer index. The next step is to apply spherical spreading 
and the characteristic impedance correction. These are derived from the 
condition of conservation of acoustic power within a conical ray tube. 

Figure 2 is a schematic diagram of a conical ray tube of solid 
angle dfl from a source to an observer. The acoustic power, which is 
the product of the intensity and the cross-sectional area, must be con- 
stant from a source radius rg to an observer radius r^. This is 
expressed in equation form as 


<P 2 <P^<^o» 2 

r_ di2 = — : r dS2 


P<=c 


P<=r 


( 1 ) 


Expressing equation (1) in dimensionless form yields 




^yp)* ^ 

Pc(y )* 

^ o 



( 2 ) 


where pc^ and pc* are determined from the atmospheric properties table 
for the observer height and source altitude, respectively. The dimension- 
less height at the source and observer are given by 


y 


o 



(3) 


and 


^s ^ sin Y + h)/RTj. ( 4 ) 

The observer distance r and observer height h are found from the 
geometry table. ^In equation (2), ^sl^o spherical spreading 

term and pcCy^) /pciy^) is the characteristic impedance correction. 
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Subband Division 


o * 

The mean-square acoustic pressure <P is expressed as a function 

of the frequency f. The table <p^(f,t,o)> has values of frequency 
corresponding to the ANSI standard 1/3-octave-band frequencies. As shown 
in the reference, the accurate application of atmospheric absorption and 
ground effects requires that the 1/3-octave bands be divided into 
N. = 2m + 1 subbands, where m is an integer greater than 0. The 
default value of m is 2, The ratio of subband center frequencies is 


- w . ( 5 ) 


j is the index number of the subband center frequencies. The 
index number j is related to the index number of the original 
1/3— octave center frequencies i by the relation 

j = (i - l)Nj^ + h (h = 1, 2, ..., Nj^) (6) 


Then, each subband center frequency is expressed in terms of a 1/3-octave- 
band frequency as 


h-m-1 

f j = w f . 


(h = 1, 2, 


i = 1, 2, . . . , k) 


(7) 


where k is the total n\imber of 1/3-octave bands. The 1/3-octave and 
subband frequencies for m = 2 are given in table II. 

The value of the mean-square acoustic pressure for each subband is 
determined from the slopes of the 1/3— octave— band spectrum. For each 
1/3-octave center frequency f^, the slope of the spectrum in the lower 
half of the band is 


Ui = 


<p^>i 

<p^>Li 


( 8 ) 


cind the slope for the upper half is 




= u 


i+1 


( 9 ) 


The term <p^>* is the value of the mean-square acoustic pressure for 
f^. The terms^ Uj^ and Vj^ are not defined by equations (8) and (9) ; 


5.1-6 



therefore, the end slopes are established by setting and 

vjc " Then the value of the mean-square pressure for each subband 

* 

^p^^j is given by 



i <p2>iAi 

(<p2>*/Ai)vih-"‘-l 

V. 


(h — 1, 2/ •••, m) 

(h = m + 1) \ (m > 0) (10) 

(h = m + 2, ...f Nj^) 


The subband adjusting factor Aj^ is defined such that 


f <p2>* - <p2>* 


(11) 


Substituting equation (10) into equation (11) and solving for yields 


^i 


_ ^ / h-m-1 ^ h\ 

1 * 2 (“i “i) 

h=l 


(12) 


Defining the subband adjusting factor in this manner ensures that 

the sum of the mean-square acoustic pressure for the subbands equals the 
total mean-square acoustic pressure of the 1/3-octave band. 


Atmospheric Effects 

The atmospheric absorption and the ground effects are now applied to 
the subband data. The atmospheric absorption coefficient y is a 
function of frequency f and altitude y as computed by the Atmospheric 
Absorption Module (ABS) . The average absorption coefficient from the 
source to the observer is defined as 


y(f,ys,yo) = ~ 


1 

- iy. ^ 


(y) dy = 


yg ~ yp 

ys - yo 


(13) 


where y^ and y^ are given by equations (3) and (4), respectively. 
Then, the mean-square pressure with the atmospheric absorption effect 
applied is 
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(14) 


^ * 

<P^X 


j ,abs 


= <p2X 


- [2y (f ,yg,yQ)f/cJ (r-rg) 
e 


The subscript cibs indicates that the absorption has been included in the 
mean- square pressure. 


Ground Effects 

Similarly, the ground effects are applied to the subband data. The 
ground effects factor G is a function of path-length difference kAr, 
cosine of the incidence angle cos 0, dimensionless frequency Hr and 
source- to- image distance ^^ 2 ^ discussed in the Ground Reflection and 
Attenuation Module (GRA) . Figure 3 shows the source- to-observer geometry 
for the ground effects computation. Referring to the figure, the 
(quantities r, h, and Y are obtained from the geometry table. Then, 
from the Law of Cosines, the source- to- image distance is given by 


+ (2h)^ - (2r) (2h) cos (90° + Y) (15) 

or 

+ 4h^ + 4rh sin y (16) 


Then the path-length difference Ar is 


Ar = r^ r 


(17) 


and the cosine of the incidence angle is 


cos 0 


r sin Y + 2h 
^2 


(18) 


The dimensionless frequency n and the wave number k are given by 


n 


2TTPaf 

O 


and 


(19) 


k 


2lTf 


( 20 ) 
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where and are evaluated from the atmospheric properties table 

at the observer altitude y^. The value of G is now computed using the 
method presented for the Ground Reflection and Attenuation Module. The 
mean-square pressure with ground effects is then given by 

<21 


The subscript gr indicates that the ground effect has been included in 
the mean-square pressure. 


Subband Combination 

The final step is to recombine the subbands into 1/3-octave bands. 
The mean-square pressure is summed over each 1/3-octave band by the 
relation 


<p^>: = 2 <p2>; 


h=l 


( 22 ) 


where j=m(i-l) +h. Once the processing of the mean- square pressure 
has been completed for all values of reception time and observer index, 
the received noise data table is complete. 

The Propagation Module always performs the spherical spreading and 
characteristic impedance change. The user has four options concerning 
the application of atmospheric attenuation and ground effects: 

1. No atmospheric attenuation or ground effects 

2. Atmospheric attenuation effects only 

3. Ground effects only 

4. Both atmospheric attenuation and ground effects 

The appropriate option is selected by the user depending on the nature of 
the problem. 

The standard output form for the noise data at the observer is the 

* 

dimensionless mean-square pressure ^p"^^ . The user may also request 
printed output of the sound pressure level SPL in decibels, defined as 
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+ 197 


(23) 


SPL = 10 log^Q <p2>* + 20 log^^ p*(c*) 

where p* and c* are determined from the atmospheric properties table 
a. a 

at the observer altitude y^. 


REFERENCE 

1. Monteganir Francis J.: Computation of Atmospheric Attenuation of 

Sound for Fractional-Octave Bands. NASA TP-1412 , 1979. 
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TABLE I.- RANGE AND DEFAULT VALUES OF INPUT PARAMETERS 


Input 

parameter 

Minimum 

Default 

Maximum 

Nb 

1 

5 

9 

^s' 

0.1 

1 

100 

a , kg/s-m^ .... 

5 X 10^ 

2.5 X 10^ 

1 X 10^ 












TABLE II.- ANSI STANDARD OCTAVE AND 1/3-OCTAVE BANDS 


WITH 1/15-OCTAVE-BAND FREQUENCIES 


Octave 

1/3 

octave 

1/15 

octave 

Octave 

1/3 

octave 

1/15 

octave 

Octave 

1/3 

octave 

1/15 

octave 



45.7 



288 



1 820 



47.9 



302 



1 910 


50 

50.0 


315 

315 


2 000 

2 000 



52.5 



331 



2 090 



55.0 



347 



2 190 



57.5 



363 



2 290 



60.3 



380 



2 400 

63 

63 

63.0 


400 

400 


2 500 

2 500 



66.1 



417 



2 630 



69.2 



437 



2 750 



72.4 



457 



2 880 



75.9 



479 



3 020 


80 

mBM 

500 

500 

500 


3 150 

3 150 



83.2 



525 



3 310 



87.1 



550 



3 470 



91.2 



575 



3 630 



95.5 



603 



3 800 





630 

630 


4 000 

4 000 



105 



661 



4 170 



110 



692 



4 370 



115 



724 



4 570 






759 



4 790 

125 

125 

125 


800 ' 

800 


5 000 

5 000 






832 



5 250 



138 



871 



5 500 



145 



912 



5 750 



151 



955 



6 030 


160 

160 

1000 

1000 

1000 


6 300 

6 300 



166 



1050 



6 610 



174 



1100 



6 920 



182 



1150 



7 240 



191 



1200 



7 590 



1 


1250 

1250 


8 000 

8 000 






1320 



8 320 



219 



1380 



8 710 



229 



1450 



9 120 



240 



1510 



9 550 

250 

250 

250 


1600 

1600 


10 000 

10 000 



263 



1660 



10 500 



275 



1740 



11 000 
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VO 
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Figure 1.- Mapping of reception time values to source time line 
demonstrating retarded time effects. 
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Figure 3.- Source-tO“Observer geometry for computation of ground effects. 



5,2 GENERAL SUPPRESSION MODULE 
INTRODUCTION 

The General Suppression Module applies a noise suppression factor to 
a noise table produced by any ANOPP noise source module. The noise sup- 
pression factor is supplied by the user as a function of frequency, polar 
directivity angle, and azimuthal directivity angle. The output table of 
suppressed noise is in the seune format as the input noise table, 

SYMBOLS 

ambient speed of sound, m/s (ft/s) 
frequency, Hz 

2 4 

mean- square acoustic pressure, re p c 

00 00 

reference pressure, 2 x lo”^ Pa (4.177 x Ib/ft^) 

suppression factor 
polar directivity angle, deg 
ambient density, kg/m^ (slugs/ft^) 
azimuthal directivity angle, deg 

INPUT 

A table of mean-square acoustic pressure and a suppression factor 
must be provided. In addition, ambient conditions are required for compu- 
tation of sound pressure levels. 

Noise liable 

f frequency , Hz 

9 polar directivity angle, deg 

0 azimuthal directivity angle, deg 

mean-square acoustic pressure, re p^c^ 

00 00 

Ambient Conditions 
ambient speed of sound, m/s (ft/s) 
ambient density, kg/m^ (slugs/ft^) 
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Suppression Factor Table 

f frequency f Hz 

0 polar directivity angle, deg 

(j) azimuthal directivity angle, deg 

S(f,9,<j)) suppression factor 


OUTPUT 

The output is the suppressed noise table. 


Suppressed Noise Table 

f frequency, Hz 

0 polar directivity angle, deg 

(J) azimuthal directivity angle, deg 

* . 2 - 
^p^ (f ,9,(J))^ suppressed mean-square acoustic pressure, re p^c 


METHOD 


The suppression factor is defined as 


S = 


<A 

<p^> 


( 1 ) 


■K 

where <p2>! is the suppressed mean— square acoustic pressure and ^p ^ 

is the Xinsuppressed mean-square acoustic pressure. The module multiplies 
each element of the input noise table by the appropriate value of the 
suppression factor to yield the suppressed noise table. In addition, 
printed output is available of the suppressed sound pressure level SPL 

defined as 


SPL = 10 log^^Q 


10 


00 

Pref 


(2) 
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6. RECEIVED NOISE 




6.1 NOISE LEVELS MODULE 


INTRODUCTION 

Previous modules within ANOPP have been concerned with Level III 
data; that is, data that have frequency, time, observer, and source as 
independent variables. In the Noise Levels Module (LEV) , these Level III 
data are integrated over frequency to produce Level II data. Level II 
noise data are a function of time, observer, and source and are referred 
to as noise levels. 

There are many varied scales that are used to describe noise levels. 
The overall sound pressure level (OASPL) is a simple integration of the 
spectrum. The A-weighted sound pressure level (L;^) and the D-weighted 
sound pressure level (Lj^) add a weighting factor to account for the rela- 
tive annoyance of the frequency bands. The perceived noise level (PNL) 
adds a weighting factor that is a function of frequency and mean-square 
acoustic pressure. The tone-corrected perceived noise level (PNLT) adds 
to PNL the impact of sounds which have pure tone characteristics. 

The Noise Levels Module can produce tables of OASPL, L;^, Lp, PNL, 
or PNLT as a function of time and observer. One, some, or all of these 
noise level scales can be selected by the user. The PNLT table is 
required by the Effective Noise Module (EFF) to compute effective per- 
ceived noise level. This module can also sum the noise data from two or 
more sources prior to computation of noise levels. 


^k'®k 

C 


c 


ADi 


F 


f 

i 


N 


SYMBOLS 

constants 

discrete frequency correction 

speed of sound, m/s (ft/s) 

second difference of band sound pressure 

tone correction difference 

frequency, Hz 

standard 1/3-octave-band niamber 
A-weighted sound pressure level, dB(A) 
D-weighted sound pressure level, dB(D) 
noisiness index, noys 
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total noys value 


OASPL overall sound pressure level, dB 

o observer index 

PNL perceived noise level, PNdB 

PNLT tone-corrected perceived noise level, PNdB 

<p2> mean-square acoustic pressure, Pa^ (Ib^/ft^) 

SPL sound pressure level, dB 

t time, s 

W weighting function 

p density, kg/m^ (slugs/ft^) 

Subscripts: 
a ambient 

r standard sea level value 

Superscript: 

* dimensionless quantity 

INPUT 

The input is a table of the mean-square acoustic pressure as pre- 
pared by the Propagation Module (PRO) • 

Received Noise Data Table 

f frequency, Hz 

t reception time, s 

o observer index 

c* (o) speed of sound at observer, re c^ 

p*(o) density at observer, re p^. 

* 2 4 

<p2(f,t,o)^ mean-square acoustic pressure, re 
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OUTPUT 


The output of this module is one or more of the following tables, 
depending on the desires of the user: 

Overall Sound Pressure Level Tcible 
t reception time, s 

o observer index 

OASPL(t,o) overall sound pressure level, dB 

A-Weighted Sound Pressure Level Table 
t reception time, s 

o observer index 

LA(t,o) A-weighted sound pressure level, dB(A) 

D-Weighted Sound Pressure Level Table 
t reception time, s 

o observer index 

LQ(t,o) D-weighted sound pressure level, dB(D) 

Perceived Noise Level Table 
t reception time, s 

o observer index 

PNL(t,o) perceived noise level, PNdB 

Tone-Corrected Perceived Noise Level Table 
t reception time, s 

o observer index 

PNLT(t,o) tone-corrected perceived noise level, PNdB 

METHOD 

The method for each noise level scale is presented as a separate 
section. Further details and a comparison of the noise level scales are 
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presented in reference 1. If two or more received noise data tables are 
input r they are summed element by element prior to the noise level 
computation. 


Overall Sound Pressure Level (OASPL) 

The overall sound pressure level is a simple integration of the fre- 
quency spectra. It is expressed in the units of decibels in the form 

n ^ * * 2 

OASPL = 10 logj^Q y 

i=l 

★ 

where n is the number of frequency terms in the input table and ^p^X 
is the dimensionless mean-square pressure for the ith frequency band. 


A-Weighted Sound Pressure Level (L^) 

The A-weighted sound pressure level applies a weighting factor to 
each 1/3-octave band that is representative of the degree of annoyance of 
each frequency. The values of the A-weighting factors are given in 
table I and plotted in figure 1. The equation for L^ is 

n 

L^ = 10 log^Q ^ 20 log^Q Pa(^a) 

where W(i,A) is the A-weighting factor. 


D-Weighted Sound Pressure Level (L^^) 

The D-weighted sound pressure level is similar to L^ except a dif- 
ferent weighting factor is applied to the data. The values of the 
D-weighting factors are given in table II and plotted in figure 2. The 
equation for L^ is 


Ld = 10 logio I [<p2>* w(i,D^ + 20 

i=l 

where W(i,D) is the D-weighting factor. 


Perceived Noise Level (PNL) 

The perceived noise level scale applies a weighting factor, which is 
a function of both frequency and sound intensity, to the mean-square 
acoustic pressure. This weighting is accomplished with the use of a 
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noisiness rating measured in noys. Figure 3 shows the values of the 
noisiness rating as a function of frequency and sound pressure level 
(SPL) . The SPL is defined as 


SPL = 10 log^Q <p2> + 20 log^Q 


p!«)' * 


(4) 


The noisiness rating has been expressed in functional form for ease of 
computation. The noisiness rating for a given value of frequency and 
sound pressure level is 


rr 


N = { 


10 


10 


10 


10 



(SPL < Aj^) 

(SPL-Aj^)-! 

(Aj^ $ SPL ^ A2> 

B2(SPL-A2) 

(A2 ^ SPL < A3) 

Bj (SPL-A^) 

(A3 5 SPL ^ A4) 

B5(SPL-A5) 

(A4 ^ SPL ^ 150) 


(5) 


The coefficients and are functions of frequency as given in 

table III. 

The computation of PNL uses values of N. The process is 

1. Determine the value of N, in noys, from equation (5) for each 

band 

2. Determine the maximum noys value N 

^ max 

3. Compute the total noys value from the equation 

24 


Nt = + 0-15 


i nI - N, 


max 


>i=l 


( 6 ) 


for the 1/3-octave bands 

4. Compute the perceived noise level in units of PNdB from 


PNL = 40 + 33.22 log^^ 


(7) 
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Tone-Corrected Perceived Noise Level (PNLT) 


The tone- corrected perceived noise level is perceived noise level 
modified for the impact of pure tone content of the noise spectra. Pure 
tones provide an additional irritation not found in broadband noise. The 
method provides for detection of pure tone content in 1/3-octave-band 
spectra and correction of the PNL for the impact of the pure tones. 

The procedure may be illustrated with reference to figure 4 , which 
shows an example spectrum, and table IV, which gives the spectrum in 
tabular form. The following steps are performed: 

1. Compute the second difference ADj^ of the SPL, which is 


Ad. = SPL.^, - 2SPL. + SPL. - 
1 1+1 1 1-1 


( 8 ) 


2. If AOj^ < -5, check to see if SPL^ is a local maximum, that is, 
if 


SPL^ > SPL^_^ and SPL^ > SPL^^j^ 


( 9 ) 


3. If SPL^ is a loca l maximum, compute the average, or background, 
noise level SPLj^ defined as 


SPL^ = 


SPLi+i + SPL._^ 


( 10 ) 


These average values are given in table IV at frequencies of 
250 and 2500 Hz. The background noise levels are shown as 
dashed lines in figure 4. 

4. The difference between the local maximum and the background noise 
level F is 


F = SPL^ - SPL^ 


( 11 ) 


5. The discrete frequency correction C(f,F) can now be determined. 
For frequencies in the range 500 < fj^ < 5000, 


0 

C(f,F) = / F/3 
6.7 


(F < 3) 
(3 ^ F < 20) 
(20 ^ F) 


( 12 ) 
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and for frequencies in the range fj^ ^ 500 or fj^ ^ 5000, 


r 

0 

C(f,F) = < F/6 


^ 3.3 


(F < 3) 

(3 < F < 20) (13) 

(20 ^ F) 


6* The tone correction is the maximum value of the discrete fre- 
quency correction example occurs at 2500 Hz 

with a value of 2. Then the tone-corrected perceived noise 
level PNLT is computed from the PNL as 


PNLT = PNL + 


(14) 


REFERENCE 

1. Edge, Philip M. , Jr.? and Cawthorn, Jimmy M. ; Selected Methods for 
Quantification of Community Exposure to Aircraft Noise. NASA 
TN D-7977, 1976. 


6.1-7 



TABLE I.- WEIGHTING FUNCTION FOR A-WEIGHTED SOUND PRESSURE LEVEL 


1/3-octave- 
band center 
frequency 

W(i,A) 

dB 

correction 

1/3-octave- 
band center 
frequency 

W(i,A) 

dB 

correction 

50 

0.00096 

-30.2 

1 

000 

1.0 

0 

63 

.0024 

-26.2 

1 

250 

1.148 

. 6 

80 

.0056 

-22.5 

1 

600 

1.259 

1.0 

100 

.0123 

-19.1 

2 

000 

1.318 

1.2 

125 

.0245 

-16.1 

2 

500 

1.349 

1,3 

160 

.0457 

-13.4 

3 

150 

1.318 

1,2 

200 

.0813 

-10.9 

4 

000 

1.259 

1.0 

250 

.138 

-8.6 

5 

000 

1.112 

,5 

315 

.219 

-6.6 

6 

300 

.977 

-.1 

400 

.331 

00 

1 

8 

000 

.776 

-1.1 

500 

.479 

-3.2 

10 

000 

,562 

-2.5 

630 

.646 

-1.9 

12 

500 

.372 

-4.3 

800 

.832 

-.8 















TABLE II.- WEIGHTING FUNCTION FOR D-WEIGHTED SOUND PRESSURE LEVEL 


1/3-octave- 
band center 
frequency 

W(i,D) 

dB 

correction 

1/3-octave- 
band center 
frequency 

W(i,D) 

dB 

correction 

50 

0.0525 

-12.8 

1 

000 

1.0 

0 

63 

.0813 

-10.9 

1 

250 

1.585 

2.0 

80 

.126 

1 

• 

o 

1 

600 

3.090 

4.9 

100 

.191 

-7.2 

2 

000 

6.166 

7.9 

125 

.282 

-5.5 

2 

500 

11.482 

10.6 

160 

.398 

-4.0 

3 

150 

14.125 

11.5 

200 

.550 

-2.6 

4 

000 

12.882 

11.1 

250 

.692 

-1.6 

5 

000 

9.120 

9.6 

315 

.832 

-.8 

6 

300 

5.754 

7.6 

400 

.912 

-.4 

8 

000 

3.548 

5.5 

500 

.933 

-.3 

10 

000 

2.188 

3.4 

630 

.891 

-.5 

12 

500 

.724 

-1.4 

800 

.871 

-.6 



















TABLE III.- CONSTANTS REQUIRED FOR COMPUTATION 


OF PERCEIVED NOISE LEVEL 


1/3-octave- 
band center 
frequency, 
Hz 

1 

®i 

^2 

^2 

^3 

^3 

^4 

®5 

^5 

50 

49 

0.079520 

55 

0.058098 

64 

0.043478 

91.01 


52 

63 

44 

.068160 

51 

.058098 

60 

.040570 

85.88 


51 

80 

39 

.068160 

46 

.052288 

56 


87.32 


49 

100 

34 

.059640 

42 

.047534 

53 

.036831 

79.85 


47 

125 

30 

.053013 

39 

.043573 

51 

.035336 

79.76 


46 

160 

27 

.053013 

36 

.043573 

48 


75.96 


45 

200 

24 

.053013 

33 

.040221 

46 


73.96 


43 

250 

21 

.053013 

30 

,037349 

44 


74.91 


42 

315 

18 

.053013 

27 

.034859 

42 


94.63 


41 

400 

16 

.053013 

25 

.034859 

40 

.030103 

100.00 


40 

500 

16 

.053013 

25 

.034859 

40 


100.00 

.030103 

40 

630 

16 

.053013 

25 

.034859 

40 


100.00 1 


40 

800 

16 

.053013 

25 

.034859 

40 

.030103 

100.00 


40 

1 000 

16 

.053013 

25 

.034859 

40 


100.00 


40 

1 250 

15 

.059640 

23 

.034859 

38 


100.00 


38 

1 600 

12 

.053013 

21 

.040221 

34 

.029960 

100.00 


34 

2 000 

9 

.053013 

18 

.037349 

32 

.029960 

100.00 


32 

2 500 

5 

.047712 

15 

.034859 

30 

.029960 

100.00 

.029960 

30 

3 150 

4 

.047712 

14 

.034859 

29 

.029960 

100.00 

.029960 

29 

4 000 

5 

.053013 

14 

.034859 

29 

.029960 

100.00 

.029960 

29 

5 000 

6 

.053013 

15 

.034859 

30 

.029960 

100.00 

.029960 

30 

6 300 

10 

.068160 

17 

.037349 

31 

.029960 

100.00 


31 

8 000 

17 

.079520 

23 

.037349 

37 

.042285 

44.29 

,029960 

34 

10 000 

21 

.0596401 

. 29 

.043573 

m 


50.72 

.029960 

37 
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TABLE IV,- EXAMPLE PROBLEM FOR DETERMINATION 


OF DISCRETE FREQUENCY CORRECTION 
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A-WEIGHTING CORRECTION, dB 



50 100 200 500 1000 2000 

FREQUENCY. Hz 


5000 10000 


Figure l.~ Decibel correction for A-weighted sound pressure level. 
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-WEIGHTING CORRECTION, dB 
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standard One-Third Octave Band Number 
125 250 500 1000 2000 4 ( 

Frequency, Hz 

Figure 4.- Example spectrum for computation of 
tone-corrected perceived noise level. 
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6.2 EFFECTIVE NOISE MODULE 


INTRODUCTION 

The Noise Levels Module (LEV) computes various noise level scales 
which are a function of observer location and time. For consideration of 
the impact of noise due to airport terminal operations, cumulative expo- 
sure to noise is an important consideration. The Effective Noise Module 
(EFF) computes time-averaged indexes of noise exposure. 

The most commonly used noise exposure index is the effective per- 
ceived noise level (EPNL) , which is incorporated into this module. There 
are a variety of other cumulative exposure indexes; however, they require 
data from multiple take-offs and landings. Multiple take-off and landing 
problems are considered in Levels I and II of ANOPP. 

This module computes a table of EPNL as a function of observer posi- 
tion, from a table of tone-corrected perceived noise level produced by 
the LEV module. 


SYMBOLS 

effective perceived noise level, EPNdB 
number of time segments 
observer index 

tone-corrected perceived noise level, PNdB 
mean- square equivalent of PNLT 
time , s 

reception time increment, s 


INPUT 

The input is a table of the tone-corrected perceived noise level 
produced by the LEV module . 


EPNL 

n 

o 

PNLT 

<pnlt^> 

t 

At 


Input Constant 

At reception time increment, s 
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Tone-Corrected Perceived Noise Level Table 


t time 9 s 

o observer index 

PNLT(t,o) tone-corrected perceived noise level, PNdB 


OUTPUT 

The output is a tcible of the effective perceived noise level and 
initial, maximum, and final PNLT values at each observer. 

Effective Perceived Noise Level Table 
o observer index 

n(o) number of time segments 

EPNL(o) effective perceived noise level, EPNdB 

PNLTj^(o) initial PNLT value, PNdB 

PNLTjj^^^(o) maximum PNLT value, PNdB 
PNLTf(o) final PNLT value, PNdB 

METHOD 

The computation of EPNL requires the time average of the PNLT. 
First, the PNLT is converted back to mean-square pressure by the 
relation 


<p„lt2> . lofM-VlO ,1, 


The EPNL can be computed in integral form as 


EPNL = 10 log I <pnlt^> dt (2) 

10 

where and tj are the initial and final times for each segment of 

the PNLT data. The segments are determined by comparing the reception 
time increment At used in the Geometry Module to the differences in 
tabulated time. If the tabulated time entries are more than At apart, 
then a new segment is started. In addition to the value of the EPNL, the 
initial, maximum, and final PNLT values are printed for each observer to 
allow assessment of the quality of the EPNL predictions. 
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7. UTILITIES 




7.1 THERMODYNAMIC UTILITIES 


INTRODUCTION 

The prediction of aircraft-engine noise requires detailed knowledge 
of the thermodynamics of the engine. The gas properties of the fuel-air 
mixture are needed to analyze the thermodynamic processes. All data on 
engine processes are stored as total variables in ANOPP. The capability 
to convert these to static variables is needed for the noise-prediction 
modules. The engine utilities perform both of these functions within 
ANOPP. This allows increased efficiency in the design of the functional 
modules . 

A brief description of each utility in ANOPP is given in the 
following : 

1. Gas-Properties Utility - Computes the ratio of specific heats, 
specific enthalpy, and specific entropy function as a function of tempera- 
ture for given values of absolute humidity, fuel-to-air ratio, and 
ambient temperature. 

2. Flow-Variables Utility - Computes the static pressure, static 
temperature, and Mach number for given values of mass flow rate, total 
pressure, total temperature, and gas properties. 

The inputs, outputs, and method for each utility are presented in 
detail herein. The ranges of the input variables are given in table I. 

All methods are extracted from references 1 to 3; however, they are docu- 
mented in any standard gas dynamics text. All thermodynamic utilities 
are available for use in any module in ANOPP. 

SYMBOLS 

A flow cross-sectional area, m'^ (ft ) 

Cp specific heat at constant pressure, m^/K-s^ (ft^/^R-s^) 

f fuel-to-air ratio 

absolute humidity, percent mole fraction 
h specific enthalpy, mVs^ (ft^/s^) 

M Mach number 

m mass, kg (slugs) 

m mass flow rate, kg/s (slugs/s) 
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2 

p pressure. Pa (Ib/ft ) 

R dry-air gas constant, m^/K-s^ (ft^/^R-s^) 

R universal gas constant, kg-mVK“S^ (slugs-f t^/^R-s^) 

R gas constant, m^/K-s^ (ft^/^R-s^) 

T temperature, K (^R) 

U velocity, m/s (ft/s) 

molecular weights of gas constituents 
mass fractions of gas constituents 
Y ratio of specific heats 

p density, kg/m^ (slugs/ft^) 

cf) entropy function, m^/K-s^ 

Siabscripts : 
a air 

r reference value (see table III) 

s static 

t total 

oo ambient 

Superscript: 

* dimensionless quantity 

GAS-PROPERTIES UTILITY 
Input 

f fuel-to-air ratio 

H absolute humidity, percent mole fraction 

a 

T* conbient temperature, re T^ 

00 

T* local temperature, re 
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Output 


The gas“properties table can be interpolated to provide the gas- 
properties data when needed within ANOPP. The gas constant is also 
provided . 

R* gas constant, re R 

Gas-Properties Table 
T* temperature , re T^ 

h*(T*) specific enthalpy, re Rt^ 

(J)*(T*) specific entropy function, re R 

Y(T*) ratio of specific heats 


Method 


The constituents of the gas within the flow are those of dry air 
(nitrogen, oxygen, and argon) and the products of combustion (carbon 
dioxide and water) . The amounts of the other gas constituents are 
assumed to be negligible* Teible II gives the five gas constituents and 
their molecular weights. The mass of each gas constituent is given by 


mi 


0.75558 

^2 


0.23154 - 3.43185f 

m 3 

= ">a 

3.13753f 

"'4 


0.00622H^ + 1.28432f 

ms 


0.01289 



— 


( 1 ) 


The amount of fuel that can be burned is limited to the availability of 
oxygen in the air* Therefore, the value of the fuel-to-air ratio f 
cannot exceed 0.06767. The constants used in equations (1) are sum- 
marized in table III. Then, the mass fraction of each constituent 

gas is 



( 2 ) 
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where the total mass m is 


5 

m = ^ (3) 

i=l 

The dimensionless gas constant R for each constituent gas can be 
computed from the universal gas constant R as 




— — " 



R/W^ 




^2 

_ 1 

R/W 2 


R 


• 


• 

• 


• 

★ 





R/W^ 


(4) 


where R is the dry-air gas constant and is the molecular weight of 

the ith gas constituent given in table II. The dimensionless gas con- 
stant of the mixture is then given by the following matrix product: 



The dimensionless gas constant R is given as a function of fuel-to-air 
ratio and absolute humidity in table IV. 

The dimensionless specific heat at constant pressure c* (which is 
computed from c^ = expressed as follows: 


* . 

c^(T 


*) = R*X 2 , R*xJ 


Cp^l(T/Tr) 



(T/Tj.) 


(6) 
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where X. and R* are the component mass fractions and gas constants 
from equations ( 2 ) and (4) and T/T^ = T*T*, The values of the component 
specific heats at constant pressure given in table V. Finally, 

the ratio of specific heats Y is 


* 



The ratio of specific heats as a function of temperature and fuel-to-air 
ratio is plotted in figure 1. 

The enthalpy per unit mass of a fluid is defined as 



In general, the specific enthalpy is a function of the fuel-to-air ratio 
and the absolute humidity. Then the specific enthalpy h^ is expressed 
in terms of the component specific enthalpies of table V as 
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The specific enthalpy is plotted as a function of temperature and fuel— to— 
air ratio for zero absolute humidity in figure 2. 


The entropy function (p per unit mass of a fluid is given by 


<t> = 



T 


dT 


(13) 


which may be expressed in dimensionless form as 



(14) 


Equation (14) can be expressed in terms of T/T^. as 


<)>*(T*) 



* 



(T/Tj.) 


d(T/Tj,) 



* 


d(T/Tj,) 


d(T/Tj.) 


(15) 


or 


P* (T*) 


K( 


T T 


) - <(t:) 


(16) 


In general, the specific entropy function is a function of the fuel—to—air 
ratio and absolute humidity. Then the specific entropy function <p* is 

expressed in terms of the component specific entropy functions of table V 
as 


*^r j^l^l' ^2^2' ■■■' ^5^^ 


* 
r,l 
* 
r,2 


"r,5 


(17) 


The specific entropy function is plotted as a function of temperature and 
fuel-to-air ratio for zero absolute humidity in figure 3. 

The gas properties table is formed from equations (7), (12), and 
(17). It can be interpolated to determine Y, h* , or ())* for a given 
value of T*. Alternately, the temperature T* can be determined for 
a given value of h* or ({)* . 
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FLOW-VARIABLES UTILITY 


Input 

The inputs to this utility are the gas-properties table 
and the total flow variables. 


• * 


t 

* 




Total Flow Variables 
total mass flow rate, re Ap^ 
total pressure, re p^ 
total temperature, re T^ 


Gas-Properties Table 
T* temperature, re T^ 

h*(T'*') specific enthalpy, re Rt^ 

(])*{T*) specific entropy function, re R 

y(T*) ratio of specific heats 


Output 

The outputs are the static flow variables. 


M 



★ 

T 

s 


Static Flow Variables 

Mach number 

static pressure, re p^ 
static temperature, re T^ 


Method 


The computation can be performed with either constant ratio of 
specific heats or variable ratio of specific heats. 


Constant Ratio of Specific Heats 

As derived in reference 2, the total mass flow rate m^ is related 
to the Mach number as follows: 
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( 18 ) 



where the ratio of specific heats y is determined from T* using the 
gas-properties table. Given the mass flow rate, the Mach n^iber cannot 
be determined in closed form from equation (18) . In addition, the equa- 
tion has two roots, one for subsonic and one for supersonic flow. The 
subsonic root is the one of interest here, so the interval-halving tech- 
nique (ref. 4) is used to find the Mach number value. The Mach number is 
plotted versus the dimensionless total mass flow rate in figure 4. 

From the energy equation for the adiabatic flow of a perfect gas 
with constant specific heats, the expression for the stagnation- 
temperature ratio (ref. 2) is 


★ 



* 

T 

s 


1 + 



(19) 


Rearranging yields 



which is plotted in figure 5 for y = 1.4. 

Similarly, for an isentropic process, the stagnation-pressure ratio 
is related to the stagnation- temperature ratio by the following relation: 


Y-1 



( 21 ) 


After substituting equation (19) and rearranging, the static pressure p* 
is 2 


P 


* 

s 







Y 

Y-1 


( 22 ) 


which is also plotted in figure 5 for y = 1.4. 


7.1-8 



Variable Ratio of Specific Heats 

Computation of the static temperature and static pressure using a 
variable ratio of specific heats requires two relations. The first is 
derived from continuity and the first law of thermodynamics, and the 
second is derived from the second law of thermodynamics. The two equa 
tions are solved simultaneously. 

The continuity relation is 


m = PgUA 

and, from the first law of thermodynamics, 

U = vj2(h^ - hg) 

Substituting equation (24) into (23) and applying the ideal-gas law 
yields 


m = ^ A \/2(h^. - hg) 
s 

By rearranging and expressing in dimensionless form, equation (25) 
becomes 



From the second law of thermodynamics, the change in entropy As 
given by 


As 



- 


s 


where is the entropy function. For an isentropic process, As - 

by rearranging yields the following: 


In 



- <»s) 

R 


( 23 ) 


(24) 


(25) 


(26) 

is 


(27) 
, so 


(28) 
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By taking the exponential of both sides and putting in dimensionless form, 
equation (28) becomes 


* 



exp [-(<{)* - 


(29) 


Equations (26) and (29) can be solved simultaneously for the static 
temperature and static pressure. Two roots exist, one for subsonic flow 
and one for supersonic flow. To ensure that the subsonic flow case is 
determined, the equations are solved using an interval-halving technique. 
The static-to-total temperature and pressure ratios computed with both 
constant and variable ratios of specific heats are given in tables VI 
and VII. 

From the continuity equation (eq. (23)), the Mach number can be 
expressed as 


M 


m 

PsA^/y^ 


(30) 


Applying the perfect-gas law and expressing in dimensionless form yields 


M 



(31) 


Table VIII compares the Mach numbers for constant and variable ratios of 
specific heats. 
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TABLE I.- RANGES OF INPUT PARAMETERS 


Input 

parameter 

Minimvim 

Maximum 

f 

0 

0.06767 

percent . . . 

0 

4 

. * 

m^ 

0 

0.6847 

* 

Pt 

0.1 

10.0 

* 

Ta 

0.8 

1.2 

* 

T 

t 

0.5 

7.0 








TABLE II.- GAS CONSTITUENTS 


Index 

Constituent 

Molecular weight 

1 

Nitrogen (N 2 ) 

28.01340 

2 

Oxygen (O 2 ) 

31.99880 

3 

Carbon dioxide (CO 2 ) 

44.00995 

4 

Water (H 2 O) 

18.01534 

5 

Argon (Ar) 

39.94800 


TABLE III.- STORED PRIMARY CONSTANTS 


Mass fraction of nitrogen in dry air 0.75558 

Mass fraction of oxygen in dry air 0.23154 

Mass fraction of argon in dry air 0.01289 

Mass of oxygen required per unit mass of 

fuel burned 3.42185 

Mass of water produced per unit mass of 

fuel burned 1.28432 

Mass of carbon dioxide produced per unit of 

fuel burned 3.13753 

Mass fraction of water per percent mole fraction 0.00622 

Standard sea level temperature, K 288.15 

Universal gas constant, kg-m /K-s 8314.32 
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TABLE IV.- VARIATION IN DIMENSIONLESS GAS CONSTANT 


f 

R* for H^, 

percent mole fraction, of - 

0 

1.2 

0.000 

1.00016 

1.00763 

.005 

1.00032 

1.00775 

.010 

1.00048 

1.00787 

.015 

1.00064 

1.00799 

.020 

1.00079 

1.00811 

.025 

1.00095 

1.00622 

.030 

1.00110 

1.00834 

.035 

1.00125 

1.00845 

.OAO 

1.00140 

1.00857 

.045 

1.00155 

1.00868 

.050 

1.00170 

1.00879 

.055 

1.00184 

1.00690 

.060 

1.00199 

1.00901 

.065 

1.00213 

1.00912 
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TABLE V.- Concluded 
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TABLE VII.- COMPARISON OF STATIC -TO- TOTAL PRESSURE RATIO, ^OR 
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Figure 1.- Ratio of specific heats. 
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Temperature Ratio, T/T^ 

Figure 2.- Specific enthalpy. 
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Figure 3.- Specific entropy function. 
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Figure 4.- Mach nximber referred to total variables for 
constant ratio of specific heats, v “ 1.4. 
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